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PREFACE 


The  Central  Technical  Unit^  Hydrology  Branch,  Engineering  Division,  Soil 
Conservation  Service  (SCS),  U.S.  Department  of  Agriculture  (USDA),  uses 
electronic  data  processing  (EDP)  to  make  volume -duration -probability  (VDP) 
analyses  of  runoff  for  selected  watersheds.    Annual  maximum  streamflow 
data  for  periods  of  1,  3,  1,  15,  50,  60,  90,  120,  I85,  and  2^h  days  are 
furnished  to  SCS  by  the  U.S.  Geological  Survey  (USGS). 

The  EDP  programs  described  in  this  technical  paper:     (l)  Transfer  USGS 
taped  data  directly  to  punched  cards  ready  for  processing,  (2)  perform 
all  the  statistical  computations  needed  to  obtain  the  maximium  annual  vol- 
ume of  flow  at  each  selected  percent  chance  for  the  10  listed  periods, 
and  (3)  plot  the  computed  results  on  log  normal  probability  charts.  The 
principal  program  (0872)  uses  the  two-parameter  gamma  distribution  to 
compute  the  0-  to  99~percent  chance  events  unless  the  gamma  shape  factor 
exceeds  51;  then  it  uses  the  log  normal  distribution.    The  program  devel- 
ops "a  wide  range  of  statistics  that  apply  to  most  cumulative  probability- 
distribution  analyses  of  natural  data  in  the  earth-science  field. 

The  programs  have  been  developed  jointly  by  the  I/fe,thematical  and  Survey 
Applications  Section,  Washington  Data  Processing  Center,  Statistical 
Reporting  Service  (SRS),  and  the  Central  Technical  Unit,  Hydrology  Branch, 
Engineering  Division,  SCS,  USDA.    They  are  written  in  FORTRAN  II  for 
processing  through  the  IBM  707^-,  IBM  ikOl,  and  EAI  Dataplotter  Model 
3300  equipment  available  in  the  Washington  Data  Processing  Center.-^ 

This  technical  paper  includes  source  programs,  tabular  data,  instmjictions 
for  preparing  input  data,  and  some  examples  of  application.    It  is  pri- 
marily a  reference  to  the  analyses  being  made  in  the  Soil  Conservation 
Service.     It  is  intended  also  as  a  guide  for  making  similar  analyses  by 
other  technicians  who  are  familiar  vrLth  electronic  data  processing. 


Trade  names  are  used  in  this  publication  solely  to  provide  specific 
information.    Mention  of  a  trade  name  does  not  constitute  an  endorsement 
by  USDA  over  other  products  not  mentioned. 
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HYDROLOGY  STUDY--A  MULTIPURPOSE  PROGRAM  FOR  SELECTED 
CUMULATIVE  PROBABILITY-DISTRIBUTION  ANALYSES 


By  W.  H.  S amnions ,  hydraulic  engineer.  Soil  Conservation  Service 


IBM  707^  PROGRAM  ABSTRACT 

Title:  Hydrology  study — a  multipurpose  program  for  selected  cumulative 
probability-distribution  analyses 

Program:  O872 

Job  No. :  920603-0010 

Date:     September  I963 

Purpose:     To  compute  various  statistics  from  hydrologic  or  other  data. 
Equipment  specifications: 

(a)  10,000  word  7070/72/7^+  series 

(b)  Seven  tape  units: 
One  program  tape 

Two  input  tapes  (units  10  and  21) 

Four  output  tapes  (units  11,  12,  20,  and  22) 

(c)  Compiled  without  floating-point  hardware 

Source  language:  FORTRAN 

Accuracy:     Single  precision  floating  point 
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PROGRAM  0872  WRIIEUP 

A.  Title:    Hydrology  study — a  multipurpose  program, for  selected  cumulative 
probability-distribution  analyses 

B.  Authors:    Elgin  G.  Fry,  Head,  Mathematical  a.nd  Survey  A_pplications 
Section,  Washington  Data  Processing  Center,  SRS,  USDA,  and  William  H. 
Sammons,  Central  Technical  Unit,  SCS,  USDA 

C.  Description:    This  program  reads  specified  control  input  data  (described 
in  later  sections)  that  give  information  to  the  computer  concerning 

the  number  of  lines  plus  the  beginning-  and  ending-line  numbers  of  any 
desired  sets  of  ICQ  lines  for  the  specified  input.    A  maximum  of  10 
sets  can  be  requested.    A  line  is  one  horizontal  line  of  input  data  on 
one  IBM  80 -column  card.     Only  79  columns  are  used;  see  data  card  format. 

A  maximum  of  100  lines  can  be  processed  for  each  station.    Each  line 
card  is  prepared  as  outlined  on  SRS -DP -102,  Multiple -Card  Layout  Form, 
page  11.    The  program  assumes  that  the  decimal  is  placed  just  before 
the  last  position  of  each  field.    Punching  of  decimals  in  the  data 
cards  overrides  the  program  and  the  data  field  will  be  read  incorrectly. 

Input  data  as  mean  discharge  for  10  selected  periods,  i.e.,  1,  3^  1, 
15,  30,  60,  90^  120,  183,  and  27^  days,  are  arranged  in  descending  order 
of  magnitude  of  the  original  data  values .    The  values  are  converted  to 
log  to  base  10  and  then  to  inches .    The  following  values  are  written 
as  output  for  each  stream  gage: 

Original  values  from  high  to  low 

Log  to  base  10  of  original  values 

Mean  discharge  converted  to  inches  of  runoff 

Additional  values  required  for  automatic  data  plotting  are  outlined 
in  E. 

D.  Principal  statistics: 

1.  Number  of  lines  per  station  (n) 

where  X  >  0;  2  ^  N  ^  100 

N  is  printed. 

2.  Coefficient  of  variation  (CV) 

Used  with  arithmetic  normal,  extreme  value,  etc.,  for  small 
values  of  CV;  assumes  no  outliers  in  the  sample  size,  N. 
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2  V 


(X^)  -  N  (SXV 


cv 


cv  is  printed. 
5.  Arithmetic  mean        (X  BAR) 


Used  with  arithmetic  normal,  extreme  value,  gamma,  etc.,  distri- 
butions; assumes  no  outliers  in  the  sample  size,  N. 

^  _  2  X 


X  BAR  is  printed. 
h.  Log  to  base  e  of  geometric  mean    (LN  G) 


Used  with  log  normal,  log  extreme  value,  gamma,  etc.,  distribu- 
tions; outliers  may  or  may  not  be  present — in  general  none  should 
be  present. 


E  log  X. 
Log^  G  =  — f-i 

LN  G  is  printed. 
5.  Geometric  mean  (G) 


Used  with  log  normal,  log  extreme  value,  gamma,  etc.,  distribu- 
tions; assumes  no  outliers  in  the  sample  size,  N. 

G  =  e(l°Se  °) 

G  is  printed. 
6.  Ratio  of  X  to  G 


Used  in  maximum-likelihood  solutions  of  log  normal,  log  extreme 
value,  gamma,  etc.,  distributions;  ass\ames  no  outliers  in  the 
sample  size,  N. 

Z  X 


RATIO  X/G  is  printed  for  X/G. 


h 

7.  Coefficient  of  variation  (maximum  likelihood)  (CV"^) 

Used  if  the  coefficient  of  variation  is  larger  than  l/j*    Used  with 
arithmetic  normal,  log  normal,  extreme  value,  log  extreme  value, 
gamma,  etc.,  distributions;  assumes  no  outliers  in  the  sample 
size,  N. 

CV*^  -  [  (x/g)^  -  1  ] 


. • .  CV*  =  V  CV*^ 

CV*  is  printed. 

8.  Gamma -shape  parameter  (y) 

Used  in  two-  and  three-parameter  gamma  distributions  and  in  items 
9  and  10. 

Y  =  Log^(X/G) 

(a)  If  0  <  Y  <  0.5772 

y  =  Y  "^(0.5000876  +  0.l6i^8852Y  -  0 .0'5hh2^kY^) 

(b)  If  0.5772  <  Y  <  17.0 

^  8.898919  +  9.05995OY  +  0.977^373Y^ 
^      Y(17. 79728  +  ll.968i+77Y  +  Y-^) 

Values  of  Y  and  GMMA  are  printed.     If  Y  >  I7.O,  the  remaining 
calculations  for  the  respective  set  are  not  made. 

9.  IM  S.D. 

Used  if  a  biased  estimate  of  the  log  standard  deviation  by  the 
maximaim-likelihood  method  is  needed;  assumes  no  outliers  in  the 
sample  size,  N. 

LN    S.D.  =  y  2  Y 
LN    S.D.  is  printed. 
10,  LN  S.D.* 

Used  if  an  unbiased  estimate  of  the  log  standard  deviation  by 
the  maximum -likelihood  method  is  needed;  assumes  no  outliers  in 
the  sample  size,  N. 


LN    S.D.*  =  y  py/JL^ 


LN    S.D.*  is  printed. 
11.  LN  S.D.** 

Used  if  an  imbiased  estimate  of  the  log  standard  deviation  by 
the  method  of  moments  is  needed^  i.e.;  log  normal  distribution; 
assumes  no  outliers  in  the  sample  size^  N.    The  solution  is 
obtained  as  in  item  2. 


LN    S.D.**  =  CV  |LN  g| 


where  Z  =  log  X. 

^e  1 


LN    S.D.**    is  printed^  but  the  definition  is  not  printed  on 
the  computer  print-out  sheets.     See  item  l8. 

12.  Sqiiare  root  of  gamma    (SQR  GAMMA.) 
Used  in  items  15,  l6,  and  17- 

SQR  GAiyiMA  =  '^~y~ 

SQP  GAMMA,  is  printed. 

13.  Scale  parameter  (BETA) 
Used  in  item  15. 

B  =  (X/y) 

BETA  is  printed.  BETA  is  calculated  in  the  original  units .  If 
the  drainage  area  is  not  zero,  BETA  is  coded  as  BETA*,  inches. 

ih.  Scale  parameter  converted  to  inches  (BETA*) 

B*  =  [     (d)     (0.037190085)     (b)  ]/M 

where  D  =  duration  in  days 

M  =  drainage  area  in  square  miles 

if  M  =  0,  B*  =  B 

B*  is  in  inches  of  depth  and  is  printed  as  BETA*. 
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15.  BETA  times  SQP  GAMMA;,  the  standard  deviation  in  the  original 
"units 

B  *  S(GA)     =  B^^  ■ 
B  ^  S(GA)  is  printed. 

16.  BETA-^  times  SQR  GAMMA,,  the  oiased  standard  deviation  in  inches  of 
depth 

B*  (S(GA)     =  (B*)V~Y~ 
B*  (S(GA)  is  printed. 

17.  U(P,I')  times  BETA*  times  SQR  GAMMA 

U(B*(S(GA)     =  (U)     (B*)  (V~r) 
where  U=  Value  of  f (l' ,P) 
P=  (y  -  1) 

If  P  is  between  two  levels  of  the  shape  parameter  (y  "~  l)^  the  P 
value  for  selected  probabilities  of  I'   (greater  than  I'  is  a 
selected  probability  level)  is  determined  by  linear  interpolation. 

See  tables  1  and  2  for  I'  =  0.0;  0.2;  1;  2;  k;  10;  20;  50;  8O;  90; 
and  99  percent  (greater  than) . 

Table  1  is  used  if  P  =  -0.95(0.05)^    (pp.  I2-I7) . 
Table  2  is  used  if  P  = -K).05(0.5)7^(l)l6i^    (pp.  I8-50). 
Table  5  supplements  table  1    (pp.  31-^l)« 

The  computer  program  uses  an  abbreviated  table  that  includes  all 
of  table  1  and  values  for  P  =  J^-. 5(0 .5)38(1)50  from  table  2.  If 
P  >  50  (or  GAMMA  >  51);  the  log  normal  is  computed.     See  items 
11  and  18. 

U(B*(s(GA)  is  printed  for  the  11  computed  points  listed. 

18.  Log-normal  solution  (calculated  only  if  GAMMA  >  51 .O) 

X.  =  [  (e(log^G  +  K^S         X)  )     (D)   (0 . O57190083)  ]/M 

e 

N^2(Z.^)  -  N(I]Z.  )^ 
^^^z^^  =    (N  -1)     (EZ. )^ 

where  M  =  drainage  area  in  square  miles 

D  =  duration  in  days  ■ 


=  reduced  variate  for  selected  probability  levels  (item  I9) 

Z.  =  log  (X.)  is  the  log  to  base  e  of  the  transformed  variate 
1  e    1  ^ 

=  log  standard  deviation  or  LN  S.D.** 

(CV  )^  is  the  unbiased  coefficient  of  variation  o'f  the  log  trans- 
formed variate  and  is  not  printed.     See  item  11.- 

19.  Other  cumulative  probability  distributions 

Reduced  variates  for  selected  distributions 
for  given  return  periods 


Return 

'''Proba- 

y 

K 

K 

value 

period 

bility 

value 

valQe 

1.01 

99 

-1.52718 

-1.6i408 

-2.3263 

1.05 

95 

-1.09719 

-1.3055 

-1.6449 

1.11 

90 

-0.83^3 

-1.1003 

-1.2816 

1.25 

80 

-0.47588 

-0.8211 

^.84l6 

2 

50 

0.36651 

-0.1643 

0.0000 

5 

20 

0.7194 

0.8416 

10 

10 

2.25037 

1.3045 

1.2816 

25 

i^- 

3.19853 

2.0438 

1.7507 

50 

2 

3 . 90194 

2.5923 

2.0537 

100 

1 

4.60015 

3.1367 

2.3263 

500 

0.2 

6.21361 

4.3947 

2.8782 

10^ 

0.001 

11.51292 

8.5265 

4.2649* 

PM 

0.00 

21.87455 

15.82576 

5.250*^ 
6.015 

*    Not  used  in  this  program. 

Probability  (percentage  greater  than). 

Seven  digits  +5  considered  for  probability;  this  is  a  graphic 
relationship. 

Other  sources  give  these  reduced  variate  values  different  mag- 
nitudes.   For  seven  9's  and  a  5,  5.0^0  to  5.327O,  5-250  is  used 
as  a  compromise.    For  nine  9's  and  a  5^  5-485  to  6.110,  6.OI5  is 
used  as  a  compromise  since  from  a  graphic  plotting  5-250  and  6.OI5 
are  consistent. 


The  return  period  (T)  is  defined  as 


where  I'  =  the  probability  (greater  than) 

For  the  extreme  value  theory^  y  value 

Y  =  -  i-a  {-  in  1) 

where  I  =  the  probability  (less  than) 

in  =  the  natural  logarithm 

The  equation  for  computed  for  selected  y  values  (y  is  a  function 
of  the  return  period) 

X.  =  u  +  (1/a)  y 

where  u  =  the  mode 

(l/a)  =  the  logarithmic  rate  of  increase  (slope) 

y  =  the  reduced  variate  for  selected  probabilities 

or  for  the  log  extreme  value  distribution 

X,  =  log  -1  t         +  (l/a)^g  y  ] 

or  in  terms  of  K  values 
g 

X.  =  X  +  K  S 
1  -    g  X 

where  X  =  the  arithmetic  mean,,  item  3 

Kg  =  the  Gumbel  reduced  variate  for  extreme  value  distribution 

^x  ^         standard  deviation  in  the  original  units ^  which  can 
be  calculated  from  item  2  or  7^  depending  on  the  size 
of  the  coefficient  of  variation 

S    =  CV  X  or  CV^  X 

and  in  terms  of  the  log  extreme  value  distribution 

X    =  log        [  X  .    +  K    S  ,  ] 
1  ia  ~    e  ^i.p-' 


ig  ~    g  xii 
to  log  base 

^x^g  ^  ^'^^^  ^'  ■^^^  °^        converted  to  log  base  10 


where  X^^  =  item  k  converted  to  log  base  10 


9 


^^^^^  ''^i'tii  ^he  normal  and  log  normal  distributions. 

For  the  normal 

X.  =  X  +  K  S 
1         —    n  X 

For  the  log  normal 


Both  Kg  and       assijme  popiUation  limitations  where  N  =  oo. 

Sample  distributions  of  K    and  K    can  be  considered  if  needed. 

§  n 

E.  Preparation  of  input  data:    Data  must  be  p^jnched  on  80-column  cards 
(only  79  columns  are  used  in  this  program)  in  the  exact  layout  speci- 
fied on  the  multiple-card  layout  form.     Columns  1  and  2  are  used  to 
identify  each  line  of  data  arranged  from  largest  to  smallest  (YR  on 
the  print -out  exhibit,  pp.  91-93).    Identification  can  range  from  00 
to  99.    There  are  10  fields  or  columns  of  input  data— 7  fields  of  8 
digits  and  5  fields  of  7  digits.    Column  80  was  not  used  in  this 
version.     It  was  reserved  for  a  binary  code  for  the  IBM-709O  computer 
program,  which  preceded  this  TBR-'JOjh  program. 

F.  Header  card  or  parameter  card:    Data  must  be  punched  on  80-column 
cards  in  the  exact  layout  specified  on  the  multiple-card  layout  form. 
Columns  1  through  10  are  for  station  identification.    Columns  11 
through  20  are  for  drainage  area  (sqLiare  miles).    The  period  of  record 
is  always  equal  to  the  total  or  set  1  and  must  be  coded  as  001  in 
columns  21  throiigh  23  to  total  N  in  columns  2k  through  26.    Sets  2 
through  10  can  be  any  sequence  of  numbers  from  001  to  N  or  parts 
thereof.    It  is  possible  to  use  all  columns  from  21  through  80  if  10 
sets  of  computations  are  desired.     If  more  than  10  sets  are  needed,  a 
new  station  is  designated. 

G.  End-of-file  card:    All  9's  must  be  punched  in  columns  1  through  30  on 
an  80-column  card  as  shown  in  the  layout  on  the  multiple -card  layout 
form. 

H.  Data  plotter  preparation:    During  the  processing  of  each  set  or  subset, 
the  values  of  CAPX(i)  (conversion  of  mean  discharge  to  inches)  are 
converted  to  log^^  and  written  on  tape  unit  20  to  be  used  as  input  to 
program  09II,  which  prepares  the  values  for  use  with  the  automatic 
data  plotter.     If  GAMMA  ^  51.0,  the  values  of  UB*  V  y    are  converted 

■to  log^o         ^^itten  immediately  after  the  CAPX(i)  values.     If  GAMMA 
is  >  5I.O,  the  log  normal  solutions  for  X.  are  converted  to  log  10  and 
written  out. 
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The  ordinate  ranges  froin  a  lower  limit  of  10      =  0.0  or  0.0001  units  to 
no  set  upper  limit.    The  abscissa  is  the  normal  reduced  variate  (K^^) 
where  a  lower  limit  of  -if. 0000  =0.0  and  +6.0000  =  0.9999  or  1.0.  Plot- 
ting of  the  observed  ordered  data  is  based  on  normal  order  statistics 
for  the  abscissa  {^^)  values.    They  can  be  plotted  in  line  or  dot  mode. 
The  11  computed  points  are  given  in  items  I'J ,  l8  and  19  of  D.  The 
ordinate  is  the  sum  of  the  log  to  base  10  of  the  observed  or  computed 
data  and  +^.000000^  which  is  then  multiplied  by  0.1.    A  scale  range  must 
be  selected  manually  for  the  data  plotted  from  the  ADP  print -out.  See 
programs  O9IO  and  09II  in  the  exhibits . 

Print -outs  of  tables  1,  2,  and  3  and  the  different  programs  follow. 
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Table  5 '--Percentage  points  of  the  chi -square  distribution-^ 

V  =  0.1(0.1)  10.0 

Y  =  0.05(0.05)  5.0 
p  =  ^.95(0.05)  h.o 

The  values  in  the  body  of  Barter's  table  2  are  the  values  of  chi  sqiiare 
that  correspond  to  these  values  of  v  (degrees  of  freedom^  reading  across, 
then  down)  for  P  =  .0001,  .0005,   .0010,   .0050,  .0100,  .0250,   .O5OO,  .1000 
(.1000),  .9000,  .9750,  .9900,  .9950,  .9990,  .9995,  and  .9999  (probability 
(less  than)  and  I'  =  100(l  —  P)).    Table  1,  Percentage  points  of  incom- 
plete gamma -function  ratio,  I'   (U,P)  was  based  on  this  table.    Table  3  (10  pp.) 
supplements  table  1  for  probabilities  not  listed  in  table  1.  Basic 
formulas  for  this  table  follow: 

=  X./p  =  uVT"  X.  =  upVT" 

u  =  2y  3  =  X/y  X.  =  (X./p)(3) 

Example  1:    Given  >=0.25  p=Y~l=-0.75 

u  =  2y  =  0.50  =  O.875289E  01 

I'  =  0.1  or  P  =  0.9990  u  =  8.75289 

Example  2:    Given  y  =  0«50        v  =  2y  =  1.0       p  =  -O.5O 

!•  =  50.0  or  P  =  0,5000 


=  uVo3o       X   =  o.i^5^956 

0.2271^-68  =  u(  0.7071067811) 
u  =  0.32168833  or 


0.321689  as  read  directly  from  table  1, 
page  1  of  5 


After  Harter,  H.  L.  More  tables  of  the  incomplete  gamma -function  ratio 
and  of  percentage  points  of  the  chi -square  distribution.    U.S.  Air  Force 
Aerospace  Res.  Lab.. Tech.  Rpt.     6k~123.     92  pp.    196^4-.     (For  sale  Office 
of  Technical  Services,  U.S.  Dept.  Commerce,  $2,25.) 
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0, 
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02 

0 
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02 
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0. 
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02 

0. 
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02 

0. 
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02 

0, 
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02 

0. 

.  1  17807E 

02 

8000 

0. 

.  129634E 

02 

0. 
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02 

0. 

. 1 32029E 

02 

0. 

.  1 33225E 

02 

0. 

.  1 34420E 

02 

9000 

0. 

.  1 5467o£ 

02 

0, 

,  1 55979E 

02 

0. 

.  1 57278E 

02 

0. 

.  158576E 

02 

0. 

.  1 59a72E 

02 

9500 

0. 

.  1 77544E 

02 

0, 

.  1 7a92VE 

02 

0. 

.'1  8031  1  E 

02 

0. 

.  181692E 

02 

0. 

.  183070E 

02 

9750 

0. 

.  1 99020E 

02 

Oi 

200477E 

02 

0. 

.201931E 

02 

0. 

.203383E 

02 

0. 

.204832E 

02 

9900 

0. 
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02 

0, 

.227494E 

02 

0. 

.229029E 

02 

0, 
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02 

0< 

.232093E 

02 

9950 

0. 
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02 

Oi 

.2471 1 9E 

02 

0. 

.24a710E 

02 

0  i 
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02 

0. 

.251 a82E 

02 

9990 

0. 

.289064E 

02 

Oi 

.290789E 

02 

0. 

.292491E 

02 

Oi 

,  294 189E 

02 

0. 

.295883E 

02 
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0< 

.307230E 

02 

Oi 

,308978E 

02 

0. 

.31 0722E 

02 

Oi 

312462t 

02 

0. 

.314  1 96E 

02 

9999 

0. 

.34831 7E 

02 

0, 

,3501 54E 

02 

Oi 

.351 987E 

02 

Oi 

3538 16E 

02 

0'. 

.355640E 

02 
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Print -out  of  FORTRM  program  No.  I519 
by 

H.  A.  Richardson 


This  program  was  used  to  compile  table  1  and  part  of  table  2,  Percentage 
points  of  the  incomplete  gamma -function  ratio,  for  the  following  P  values: 

p  :=  -0.95  (0.05)^ 

P  =  k.^  (0.5)58 

P  =  39  (1.0)50 
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COMPILE        RUN  FORTRAN 
C  PERCENTAGE   POINTS   OF    THE    INCOMPLETE   GAMMA-FUNCTION   RATIO  FOR 

C  THE   FOLLOWING  P-VALUES. 

C  P    =  -.95(.05)4 

C  P=4.5(.5)38 
C  P    =    39 ( 1  .0  )50 

C  JOB   NO.    920603-0010  PROGRAM    NO.  1319 

C  PROGRAMMED    BY    H.A.RICHARDSON  JULY    13,  1964 

DIMENSION    A(429,5),  8(195,12) 

REWIND  1 

REWIND  2 

REWIND  3 

READ    INPUT    TAPE    1,    1»     ((A(I,J),    J    =    1,    5),     I    =    If  429) 

1  FORMAT { 1  OX ,5E 14 .6  ) 

N    =  0 

DO   2    I    =    1 ,  39 

DO   2   J    =    1  ,  11 

K    =    13   -  J 

M   =   5*  I 

L    =    M    -  4 

N    =    N    +  1 

IB    =  1 

DO    2    lA    =    L,  M 

B  (  I  A  ,  K  )    =    A  (  N  ,  I  B  ) 

2  IB    =    IB    +  1 

B ( 1 , 1 )  =  -  0.95 
DO   3    I    =2,  100 

3  B(I,1)    =   B(I    -    1,1)    +  0.05 
B( 101 , 1 )    =  4.5 
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DO    4    I    =    102.  192 
4    B(I»1)    -    B(I    -  +  0.5 

WRITE    OUTPUT    TAPE    2,  5 
50F0RMAT (IHl , 13X,95HPERCENTAGE    POINTS    OF    THE    INCOMPLETE    GAMMA-FUNCT I 
ION    RATIO    FOR    THE    P-VALUES    INDICATED    IN   COLUMN    1 /// 1  HO . 1  OX , 1  HP , 6X , 
25H1 .000,4X,5H0.998,4X,5H0,990,4X,5H0.980,4X,5H0.960,4X,5H0.900,4X, 
35H0.80  0,4X.5H0.500.4X.5H0.2  0C,4X,5H0. 100.4X ,5H0.0  1 0 ) 
WRITE    TAPE    3.     ((B(I,J).    J    =    1«     12).     I    =    1.  168) 
WRITE    OUTPUT    TAPE    2,    6.     ((B(I.J).    J    =    1.     12),     I    =    1.  168) 

6  FORMAT ( 1 HO/1 H0,4X , 1 2F9. 3/ (5X , 1 2F9. 3 ) ) 

WRITE    TAPE    3,     ((B(I.J),    J    =    li     12),     I    =    170,     192,  2) 

WRITE   OUTPUT    TAPE   2,    7,     ((B(I,J),    J    =    1,    12),    I    =    170,    192,  2) 

7  FORMAT ( 5X ,  1 2F9. 3  ) 
END   FILE  2 
REWIND  1 

REWIND  2 
REWIND  3 
TYPE  ICQ 

100  FORMAT (24HEND    OF    JOB,    UNLOAD  TAPES) 
TYPE  101 

101  FORMAT (55H0UTPUT    ON    TAPE    11    MUST   BE    PRINTED   UNDER   FORTRAN  CONTROL) 
STOP  3451 

END 

/ 

■  •  / 


/ 


Print -out  of  FORTRM  program  No.  O872 
by 

Elgin  G.  Fry 


This  print-out  is  for  the  main  program,  which  has  been  titled  Hydrology 
study— a  multipurpose  program  for  selected  cimilative  probability- 
distribution  analyses. 
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COMPILE       RUN  FORTRAN 

C  JOB  NO  9206030010     PROG  0872 

C  HYDROLOGY   STUDY  SCS 

C  JOB  NUMBER  920603-0010 

C  PROGRAM   NUMBER  0872 

C  WRITTEN  BY  ELGIN  G.    FRY  DP-SRS-USDA 

C  INPUT   TAPE    IS   10  AND  20 

C  OUTPUT  TAPE   IS   11,    12,.  AND  20 

C  REELS  ON   11   AND  12  ARE  FILLED   IN  SEQUENCE 

C  REEL   ON  20    IS    INPUT   TO  PROGRAM  0911 

C  WORK  TAPE  ON  22 

C 

C  END  OF   JOB  CARD— NINES  PUNCHED  IN  COLUMNS  1-30 


C 


OD I  MENS  I  ON   X(  100, 10) , Yl 100, 10) , I YR( 100) , XLOG( 100, 10) , C APX (  100, 10) , 
lIYEARt  100,10)  ,ISTART(  10),    lENO(lO),    DAY(IO),    N(IO),  OBS(IO),. 
2SUMXt  10)  ,.  SUMXXI  10)  ,    CVLIO),    XM(10)i    GM(IO),    SD(IO),  YX(IO), 
3BETA(10),   BETA2(ia),   GAMMAIIO),   SQGAM(IO),   UdO.ll),   P  (  10 )  ,  CZ  (  10 )  , 
4TABLE( 180, 12) ,    XLN(IO),    PK(ll),    XXN(ll),  CHANCE(ll) 

EQUIVALENCE  (X,CAPX) 

REWIND  1 

REWIND  2 

REWIND  3 

REWIND  7 

REWIND  8 

L   =  2 

NUM  =  0 

CHANCE(l)  =  0.0 
CHANCE{2)   =  0.2 


CHANCEO)    =  1.0 
CHANCEt-V).  =  2.0 
CHANCB(5)    =  A.O 
CHANCE(6)    =  10.0 
CHANCE(7)    =  20.0 
CHANCE(8)    =  50.0 
CHANCE(9)    =  80.0 
CHANCE  (  10)    =  90.0' 
CHANCE(ll)    =  99.0 
DAY(l)    =  1.0 
0AY(2)    =  3.0 
DAY(3)    =  7.0 
DAY(^)    =  15.0 
DAY(5)    =  30.0 
DAY(6)    =  60.0 
DAY(7)    =  90.0 
0AY(8)    =  120.0 
DAY(9)    =  1?3.0 
DAY(  10)    =  27^^.0 
PK(1)    =  S.-'.Sj 
PK(2)    =  2.878 
Pr( 3)    =  2.326 
PK(4)    =  2.054 
PK(5)    =  1.751 
PK(6)    =  1.282 
PK(7)    =  0.842 
PK( 8)    =  0. 0 
PK(9)    =  -0.842 
PK( 10)    =  -1.282 
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PK( 11 )    =  -2.326 

C  READ   IN   TABLE   OF   U  VALUES   FOR  VARIOUS   LEVELS  OF  PROBABILITY 

READ  TAPE   8,    ( ( T ABLE ( I , J ) ,    J=l,12),  1=1,168) 
READ   TAPE   8,    ((TABLEIUJ),    J=l,  12),    1=169,  180) 
C  READ   lU  PARAMETER,  CARD 

20READ   INPUT   TAPE    1,    ^,    6TAT1,    STAT2,.  AREA,    (ISTART(I),  lENDd), 
11=1,10) 

4  FORMAT(2A5,F10. 3,2013) 

IF    (  ISTART  (  1  )-99),  5,  IT'h,  174 
C  READ   IN  DATA  CARDS 

5  K  =  lEND(l) 
00  6    1  =  1,. K 

READ   INPUT   TAPE   1,   8,.  lYRlI),    ( X  {  UJ  ) ,    J  =  l,10) 

6  WRITE   TAPE   9,  lYR(I),    (X(I,J),  J=l,10) 

8  FORMATt I2,7F8.1,3F7.1) 
REWIND  9 

DO   170  KK=1,10 
DO   9  I=1,K 

9  READ   TAPE   9,  lYRII),    < X ( UJ ) ,  J=l,10) 
REWIND  9 

Kl  ISTART(KK) 
K2   =  lEND(KK) 
K3   =  K2  -  1 
IF    (Kl)  2,2,10 
C  CHANGE   STORAGE  LOCATIONS 

10  DO   12  I=K1,K2 
DO   12  J=l,10 
XL0GII,J)    =  0.0 
YM.,  J)   =  X(  I  ,  J) 


12  IYEAR(I.,J)    =  IY.<(I) 

PLACE   Y   VALUES   ORDERED  HIGH   TO  LOW 
DO    16   J=l , 10 

13  DO    16  I=Kl,K3 

IF    (Y[  I  ,  J)    -    Y(  in,  J)  )  l^,16^16 
1^   HOLD   =    Y( I , J) 

Y( I, J)    =   Y( I+l , J) 
Y  (  I.+  l ,  J  )    =  HOLD 
IHOLD   =    I  YEAR ( I , J) 
IYEAR(I,J)    =    I  YEAR(  Ifl,  J). 
IYEAR(  I  +  I, J)  =IHDLD 
GO   TO  13 

16  CONTINUE 

CALCULATE   LOG   TO  BASE   10  OF   Y  VALUES 

DO   18  I=Kl,K2 

DO    18  J=l,10 

IF    (Y( I, JJ )  18,13,17 

17  XLOG(  I  ,  J)    =   LOGXF  (  Y(  I  ,.J)  ) 

18  CONTINUE 

CALCULATE   CAP(X)  VALUES 
IF    (AREA)  19,19,22- 

19  DO   20  Ii=Kl,K2 
DO   20  J=l,10 

20  CAPX( I , J )    =   Y (  I i  J ) 
GO   TO  26 

22   CC   =   0.03719008   /  AREA 

DO  2^  I*K1,K2 

DO   2A  J=l,10 
24  CAPX;i,J)    =   Y(I,J)    •   DAY(J)    •  CC 
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C  WRITE   HEADER  LINES 

26   WRITE   OUTPUT   TAPE   L,    28,    STATl,    STAT2,    AREA,  KK 

280F0KMAT ( IHl , lOX , 61KHIGHEST   MEAN   DISCHARGE   GAMMA   FUNCTION  PARAMETERS 
I   FOR   STATION    ,  2  A5  ,AX  r^SHAREA   =,  ,.F  I  0 .  3 ,  5X  ,  3HS  E  T  ,  I  3  ) 
WRITE   OUTPUT   TAPE   L,  30 
300F0RMAT (  IHO, 1I5HYR  X(l)        LOG   XCI)        CAP    X(I)  YR  X( 

12)        LOG   X(2)        CAP   X(2)  YR  X(3)        LOG  X(3)        CAP  X(3)  ) 

C  WRITE   DATA  VALUES 

DO   32  I=Kl,K2 

320WRITE   OUTPUT   TAPE   L,    3'f,    IYEAR(I,.l),    Yd, I),    XLOGU,!),  CAPXd.l), 
IIYEAR(I,2),    Yd, 2),    XLnGd,2),    CAPX(I,2),    IYEARd,3),    Yd, 3), 
2XLOGd.,3),  CAPX(I,.3) 
34  FORMAT { I3,F12.4,F10.6,F10.4,2{I8,FI2.4,F10.6,F10.4)) 

WRITE   OUTPUT   TAPE   L,  36 
360F0RMAT(  IHO,  115HYR-  Xl^f)        LOG  X(4)        CAP   X{4)  YR  XI 

15)        LOG   X(5)        CAP   X(5)  YR  X(6)        LOG  X(6)        CAP   XI6)  ) 

C  WRITE   DATA  VALUES 

DO   38  I=K1,K2 

380WRITE   OUTPUT   TAPE   L,    34,    IYEARd,4),    Y(I,4),    XL0G(I,4),  CAPXd,4), 
lIYEARd,5),    Yd, 5),    XL0Gd,5),    CAPX(I,.5),    IYEARd,6),    Yd, 6), 
2XL0G(I,6),  CAPX(I,,6) 
WRITE   OUTPUT   TAPE   L,  40 
400FORMATdHO,  II5HYII  X(7)        LOG  X(7)        CAP   X(7)  YR  X( 

18)        LOG   X(8)        CAP   XIB)  YR  X(9)        LOG  X(9)        CAP   X(9)  ) 

C  WRITE   DATA  VALUES 

DO  42  I=Kl,K2 

420WRITE   OUTPUT   TAPE   L,    34,    IYEAR(I,7),    Y(I,7),    XL0G(I,7),  CAPXd,7), 
lIYEARd,8),    Yd, 8),    XL0S(I.,8),    CAPX(I,8),    IYEARd,9),    Yd, 9), 
2XL0Gd,9),  CAPX(I.,9) 
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WRITE   OUTPUT   TAPE  L, 
^4   FORMAT( 1H0,35HYR  XllO)      LOG   XllO)    CAP    X(10)  ) 

DO   46  I=K1,K2 

460WRITE   OUTPUT   TAPE   L,    34,    IYEAR(I,10),    Y{I,10),    XL0G(I,10),  CAPX(I, 
110) 

DO  41  I=K1,K2 

DO  41  J=l,10 

IF    ICAPXIUJ))    55,  55,  43 
55   CAPX( I  ,  J)    =  -9.0 

GO   TO  41 

43   CAPXd.J)    =   LOGXFICAPXl  >,  J)  ) 
41  CONTINUE 

KOOE   =   K2   -   Kl    ♦  1 

WRITE   OUTPUT   TAPE   7,    45,    STATl,    STAT2,    AREA,    KK,  KODE 
45  FORMAT  (  2A5,.F10.  3,  I  3,  I  3  ) 
DO  49  I=K1,K2 

WRITE   OUTPUT   TAPE'  7,    47,    (CAPX(Ir-J)»  J=l,10) 

47  FORMAT! 10F10.4) 

49  CONTINUE 

CALCULATE   VALUE   DF   N  FOR  EACH  DAY 
DO  48  1=1,10 

48  Mt  I)    =  0 

DO   52  I=K1,.K2 

DO   52  J=l,10 

IF    (Y(I,J))    52, 52-,  50 

50  N(J)    =   N(J)    +  1 
52  CONTINUE 

WRITE   OUTPUT   TAPE   L,    53,  AREA 
530F0RMAT( 1H1,30X,65KVOLUME-DURATION-PROBABILITY   ANALYSIS   FOR  SELECTE 
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ID   WATERSHEDS  /  L  HO ,  1  <VHN  AME   OF   S  TRE  AM ,        ,  1  3HGAGE  LOCATION/ 

21H0,8HUSGS  NO.  ,  I  5X.,  7HC  TU  NO  .  ,.1  1  X  ,.1  6HPER  I OD  OF   RECORD, 32X, 
315H0RAINAGE   AREA  =,F8.3,6H   SQ  MI    / IHO , 20HDUR AT  I  ON    IN  DAYS  r.lGX, 
^IH3,10X,IH7,9X,2H15,9X,2H30,9X,2H6  0,9X,2H90,8X,3H120,8X,3HI83,8X, 
53H274  ) 

DO   5^    I.=  l,lO 
54  OBS( I )    =   N( I) 

WRITE   OUTPUT   TAPE   L,    56^    (Nil),  I=.I,10) 
56  FORMATC  lOHON  ,  lOHl) 

CALCULATE   MEAN,  CV 

DO  58  1=1,10 

XMII)    =  0,0 

CV(  I  )    =  0.0 

CZ(  I)    =  0.0 

SD(  I  )   =  0.0 

XLNI I )    =  0.0 

GM( I )   =  0.0 

SUMXI I )    =  0.0 
58   SUMXXt I )    =  0.0 

DO  60  I=Kl,K2 

DO  60  J=l,10 

SUMXIJ)    =   SUMX(J).  ♦  Y(f,J) 

60  SUMXX(J)    =   SUMXX(J)    +   YtI,J)    «  Y(I,.J) 
DO  62  J=I,10 

IF   (OBS(J))  62,62,61 

61  XM(J)    =^  SUMXIJ)    /  OBS(J) 

OCV(J)  =  ((OBS{J)  •  OBSIJ)  ♦  SUMXX(J)  -  OBS(J)  ♦  SUMXIJ)  ♦  SUMX(J)) 
1/    KOBSIJ)    -    1.0)    •   SUMXIJ)    »   SUMXIJ)))  0.5 

62  CONTINUE 
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C  WRITE   MEANSt    CV,    AMD  SD 

WRITE   OUTPUT   TAPE   L,    6<n    (CVd),    1  =  1,10) 
6^   FORMATdOH   CV  ,10F11.5) 

WRITE   OUTPUT   TAPE   L,    66,    (XM(I),  1=1,10) 
66   FORMATdOH    X   BAR  ,10F11.4) 
C  CONVERT  LOG  AREA   TO  LOG  OF   X   TO  BASE  E 

DO  70  I=Kl,K2 
DO   70  J=1,10 
IF    (Yd, J))  70,70,69 

69  XLOG( I , J)    =   LOGF(Y( I , J)  ) 

70  CONTINUE 

C  CALCULATE   GEOMETftl.C   MEAN  AND  C{Z) 

DO  72  1=1,10 
SUMXX( I )    =  0.0 
72   SUMXI I )    =  0.0 
DO   74  I=K1,K2 
DO  7<f  J=l,10 

SUMXX(J)    =   SUMXXIJ)        XLOGd,.J)    •  XLOG{I,J) 
74  SUMX(J)    =   SUMX(J)    +  XLOGd.J) 

DO   76  J=l,10 

IF    (OBS(J).)  76,76,75 
750CZ(J)    =    ((OBS{J)    •   OBS(J)    *   SUMXX(J)    -OBSIJ)    •   SUMX(J)  »SUMX(J)) 
1/    ((OBSIJ)    -   1.0)    •   SUMX(J)    •   SUMXCJ)))    •»  0.5 

S0(J)    =   SUMX(J)    /  OBSIJ) 

XLN(J)    =  SD{J) 

CZ(J)    =   CZ{  J)    ♦   ABSFIXLNt J) ) 
GM( J)    =   EXPFCSDI J) ) 
76  CONTINUE 
C  WRI.TE   LOG  GM  AND   THE  GM 
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WRITE   OUTPUT   TAPE   L,    78,    (SD(I),  1=1,10) 
78  FORMATIIOH   LN   G  ,10FI1.6) 

WRITE   OUTPUT   TAPE   L,    80,    (GMU),  1=1,10) 

80  F0RMAT(10H  G  ,10F11.4) 

C  CALCULATE   AND  WRITE   RATIO  OF   X   BAR  /  G 

DO   82  I?I,10 
SD( I )    =  0.0 
IF   (OBS<I.))    82,82,  31 

81  SOI  I )    =   XM(  I  )    /   3M(  I  ) 

82  CONTINUE 

WRITE   OUTPUT   TAPE   L,    8^*,    (SD(I),    1  =  1,10) 

84  FORMATdOH  RATIO  X/G,10F11.5) 
C  CALCULATE   AND   WRITE   CV»   AND  Y 

DO   86  1=1,10 

YX ( I )    =  0.0 

SUMX( I )    =  0.0 

IF    (OBSJD)  86,85,85 

85  SUMX(I)    =   (SD(I)   •   SDH)   -   1,0)    »♦  0.5 
YK( I )    =   LOGF( S0( I) ) 

86  CONTINUE 

WRITE   OUTPUT   TAPE   L*    88,    <SUMX(I),.  1=1,10) 
88  FORMATdOH  CV»  ^lOFll.5) 

WRITE   OUTPUT   TAPE   L,    90,    lYX(I),  1=1,10) 
90  FORMATdOH   Y  ,10F11.6) 
C  CALCULATE  LN  S.D.   AND  LN  S.D.»   (USE   BETA  AREA) 

DO  240  I=1V10 
BETAd)   =  0.0 
IF    (OBS(n)  240,2-40,241 
241  BETAJI)    =   (YXd)   •  2.0)   »•  0.5 
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240  COMTINUE 

WRITE   OUTPUT   TAPE   L,    Z(.2,    (  BET  A  (  I  )  ,  I  =l ,  lo  ) 

242  FORMATdOH  LH   S.D.  ,10F11.6) 
DO  244  1=1,10 

IF    (OBS(I))  244,244,243 

243  BETA(I)    =   BETA(I)    .    (OBSd)    /    (OBSd)   -   1.0))  0.5 

244  CONTINUE 

WRITE   OUTPUT   TAPE   L,    246,    (BETACI),  1=1,10) 

246  FORMATdOH   LN   S.D.»  ,10F11,6) 

WRITE   OUTPUT   TAPE   L,    245.    (CZd),    I.=  l,10) 

245  FORMATdOH   LN   S.  0.      ,  1  OF  1  1  .  6  ) 

C  CALCULATE   BETA  AMD  GAMMA  AND   SOR  ROOT   OF  3AMMA 

DO   99  1=1,10 
IF    (OBSd))  247,247,248 

247  GAMMA( 1)    =  0.0 
BETA( I )    =  0.0 
SQGAMI I )    =  0.0 
GO   TO  99 

248  IF    (YXd  )-0.  5772).  92,94,94 

920GAMMAd)    =    d.o   /   YX(I)>   .    (0.5000876  +  0.1648852   •  YX(I) 
1-   0.0544274   •    YX(I)    •  YX(I)) 
GO   TO  98 
94    IF    (YXd)    -    17.0)    96,96,  170 

960GAMMA(I)    ^    (8.898919  >   9.O59950   •  YX  (I  )    +  0.9775373   .   YXd)  . 

lYXd))    /    (YXd)    .    di. 79728   *    1  1.968477   .   YX(I)    +   YX(I)    .  YXd))) 

98  BETA(  I  )    =   XMd  )    /  GAMMAd  ) 
SQGAM(I)    =   GAMMA(I)  0,5 

99  CONTINUE 
GO  TO  104 
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C  WRITE   GAMMA,    SQR  ROOT   OF   GAMMA   AND  BETA 

,104  WRITE   OUTPUT   TAPE   L,    106,    (GAMMA(I),    I  =  UIO) 
106   FORMATdOH   GAMMA  ,I0F11.5) 

WRITE   OUTPUT   TAPE   L,    108,    (SQGAM(I),  I=U10) 
108   FORMATCIOH   SQR   G AMMA , 1  OF  1 1 . 5 ) 

WRITE   OUTPUT   TAPE   L»    110,    (BETA(I),  1=1,10) 
110  FORMATdOH   BETA  ,10F11.5) 
C  CALCULATE   AND   WRITE   BETA   STAR  (PRIME) 

DO   116  1=1,10 

IF    (AREA)    114, 111^,112 

112  !F    (OBS( I ))  114,114,113 

113  BETA2(I)    =   DAYd).  »   0.03719   •   BETAd)    /  AREA 
GO   TO   1 16 

114  BETA2(  I )    =   BETA{ I) 
116  CONTINUE 

WRITE   OUTPUT   TAPE   L,    118,    (BETA2(I),  1=1,10) 
118   F0RMAT(10H  BETA*  ,10F11.5) 
C  SUMX  =   BETA   X   SQR   ROOT  GAMMA 

C  SD  =   BETA*   X   SQR  ROOT  GAMMA 

DO   120    I  =  1,T0 

SUMXd)    =   BETAd)    *  SQGAMd) 
120  SD(I)    =   BETA2(I)    •  SQGAMd) 

WRITE   OUTPUT   TAPE   L,    122,    (SUMXd),    1  =  1,10) 
122   FORMAT(10H   B   •   S ( G A ) , 1  OF  1 1 . 4 ) 

WRITE   OUTPUT   TAPE   L,    124,    CSDCI),  1=1,10) 
124  F0RMAT(10H   B*    ( S ( GA ) , 1  OF  1 1 . 6 ) 

DO    126    1  =  1, .10 
126   P( I )    =   GAMMA (  I  )    -  1.0 
C  CALCULATE   AND   STORE   U   •   B   •  S(GA) 
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DO  127  1=1,-10 
DO   127   J  =  l,.ll 

127  U( I , J)    =  0.0 
DO   152   J  =  l,.10 

IF    (OBS(J})  152,152,129 

129  IF    (P(J)-50.0)  128,128,152 

128  DO    130  1=2,180 

IF    (  P(  J)-TABLE  (  I  ,.n  )  131,1^^,130 

130  CONTINUE 
GO   TO  15? 

131  IF    (  ( P ( J)-TABLE ( I-l , 1)  )    - ( T ABLE ( I , 1 ) -P ( J )  ) )      132,  132,  133 

132  PROS  =    {P(J)    -   TABLE  (  I-l  ,  1 )  )    /    (TABLE(I,1)    -  TABLE(I-l,in 
SW  =  0.0 

GO   TO  134 

133  PROB  =    ITABLEII  ,1)-P(  J))    /    (TABLEd.l)    -  T  ABLE  (  I- 1 ,  1  )  ) 
SW  =  1.0 

134  DO    142  M=U11 

IF    (TABLE(  I-l,M*-l)    -   T  ABLE  (  I  ,  M+ 1  )  )    138,  136,  136 

136  IF    (SW)  137,135,137 

137  PROBT=   TABLE(  I  ,M«-;i  )    +    UT  ABLE  I  I  -  UM*- 1 )    -   T  ABLE  (  I  ,  M+ 1 )  )    ♦  PROB) 
GO   TO  140 

135  PROBT   =   TABLE( I-1,M*1)    -    (  I  TABLE ( I-l ,M+ 1 ) -T ABLE ( I ,  M  +  l )  )    •  PROB) 
GO   TO  140 

138  IF    (SW)  139,139,141 

139  PROBT=   TABLE( I-1,M*1)    ♦    I  I T ABLE ( I , M+ 1 )    -  T ABLE (  I - 1 , 1 ) )    •  PROB) 
GO   TO  140 

141  PftOBT  =   TA&LE(I,M*-1)    -    t  (  T  ABLE  (  I  ,      1  )   -  T  ABLE  (  I  - 1 ,  M  + 1 )  )    •  PROB) 

140  U(J^M)    =   PROBT*    SO  (  J  ) 

142  COMTINUE 
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GO   TO  152 
1^4  DO    1^6  M=l,ll 

146  U(J,M)    =   TABLE(I,M+1)    •  SO(J) 
152  CONTINUE 

WRITE   U  »   B   •    S13A)  VALUES 
WRITE   OUTPUT   TAPE   L»  I'^S 
148   FORMAT  (  1H0,.14HPE^CENT   CHANCE  |.43X,.9HU  (  B»S  (  GA  )  ) 

WRITE   OUTPUT   TAPE   Lt    154,    (CHANCE(I),    (UIJ.I),    J=l,10),  1=1,11) 
154  FORMATdH   ,.F9.  1 ,  lOF  1 1 .  4) 
WRITE   OUTPUT   TAPE   L,  161 
16I0FGRMAT ( 120H0N   =   STATION  YEARS,    CV  =   COEFFICIENT   OF  VARIATION   BY  ME 
ITHOD  OF   MOMENTS,    X   BAR   =«   ARITHMETIC   MEAN    IN   SECOND-FEET   PER  DAY,  / 
2116H  LN  G   =   NATURAL   LOGARITHUM  OF   THE  GEOMETRIC   MEANIMEAN  OF   THE  N 
3ATURAL   LOGARITHMS  OF   THE   ORIGINAL   DATA),    G  =  GEOMETRIC/    120H  MEAN 
4IN   SECOND-FEET   PER  DAY,   RATIO  X/G  =   RATIO  OF   ARITHMETIC   MEAN   TO  GE 
50METRIC   MEAN,    CV»   =  COEFFICIENT   OF   VARIATION  BY   /120H   METHOD  OF  MA 
6XIMUM   LIKELIHOOD,    Y  =   NATURAL   LOGARITHM  OF   THE   RATIO  OF   THE  ARITHM 
7ETIC   MEAN  TO  THE   GEOMETRIC   MEAN,    LN   S.D.    /114H   AND  LN   S.D.»   =  STAN 
8DARD  DEVIATION   OF   THE   NATURAL   LOGARITHMS   FOR  MAXIMUM   LIKELIHOOD,  B 
9IASED   AND  UNBIASED  RESPECT-  ) 
WRITE  OUTPUT  TAPE'  L,  185 
1850F0RMAT{ 120H   IVELY,    GAMMA  =,  SHAPE   STATISTIC,    SQR  GAMMA  =   SQUARE  ROO 
IT   OF   GAMMA,    BETA  A^NO  BETA*   =   RATIO  OF  ARITHMETIC  MEAN   TO  GAMMA  IN/ 
2116H   SECOND-FEET   PER   DAY  AND  VOLUME    IN    INCHES   RESPECTIVELY,    B   •  S( 
3GA)    AND  B»    (S(GA).  =  BETA   TIMES   THE   SQUARE   ROOT   OF  /120H  GAMMA 

41STANDARD  DEVIATION)    IN  SECOND-FEET   PER   DAY   AND  VOLUME    IN    INCHES  R 
5ESPECTI.VELY,.  U{B»S(GA)    =   VOLUME    IN   INCHES   FOR       /   51H  SELECTED  DUR 
6ATI0NS   AND  PERCENT  CHANCE   AS    INDICATED   /    IHl  ) 
DO   155  1=1,10 
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DO    155   J=l  , 1 1 

IF    (U( I , J) )    151, 151 , 153 
151   U(  UJ  )    =  -9.0 

GO   TO  155 
153  U(  I  ,.J)    =   LOGXF(U(  I  ,.J)  ) 

155  CONTINUE 

DO    166  1=1,10 

IF    (OBS( I ) )    156, 156, 157 

157  IF    ( GAMMA ( I )-51 . 0)  156,156,158 

156  WRITE   OUTPUT   TAPE    7,    164,    I,    (  U  (  I  ,.  J )  ,  J=l,in 
GO   TO  166 

C  CALCULATE    LOG-NORMAL  SOLUTION 

158  DO   160   J  =  l  ,  1 1 

XXN(J)    =    EXPFIXLNd)    +    PKIJ)    »   A&SF^  CZ  (  1)  )  ) 

IF    (AREA)  160,160,183 
183   XXN(J)    =   XXN(J)    •   DAY(I)    ♦   0.03719   /  AREA 
160  CONTINUE 

WRITE   OUTPUT   TAPE   L,    162,    I,    {XXNIJ)^  J=l,ll) 

162  FORMAT( IHO, I2,4X,11F10.6) 
DO   163   J  =  l  ,  1 1 

IF    (XXN(J))  171,171,165 
171    XXN( J)    =  -9.0 
GO   TO  163 

165  XXN<  J)    =   LOGXF ( XXN( J)  ) 

163  CONTINUE 

WRITE   OUTPUT   TAPE   7,    164,    I,    (XXN(J),  J=l,ll) 

164  FORMAT! I3,4X, 11F10.4) 

166  CONTINUE 

170  NUM  =   NUM   +    (4  »   K2   +.  37) 
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IF    (NUM-2000O)  2,172,.172 
172  END  FILE  L 

REWIND  L 

L   =  L   +  I 

NUM  =  0 

GO   TO  2 
174  END  FILE  L 

REWIND  1 

REWIND  L 

KODE   =  999 

WRITE   OUTPUT   TAPE   7,    45,'   STATl  ,    STAT2,    AREA,    KK,  KODE 
END  FILE  7 
REWIND  7 
REWIND  8 

TYPE   176  ^ 
176  FORMAT{ lOHEND  OF  JOB) 
STOP  9999 
END 


Print-outs  of  FORTRM  programs  No.  O9IO  and  No.  09II 

by 

Elgin  G.  Fry 


In  order  to  plot  on  the  automatic  data  plotter  the  computed  output  from 
program  No.  O872  must  "be  coded  and  nonnal  order  statistics  for  data- 
plotter prepared  (program  No.  09II). 

Program  No.  0910^  Conversion  of  normal  order  statistics^  was  used  in 
plotting  the  coded  data.  See  the  coded  print-outs  for  Catawba  Creek 
and  Smiths  Fork  forecast pages  91-93  and  98-IOI. 
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COMPILE        RUN      FORTRAN  1  of  1 

C  CONVERblON   OF   NORMAL  ORDER  STATISTICS 

C  JOd   NUMBER  920603-0010 

C  PROGRAM  NUM&ER   09 lO 

C  WRITTEN  BY   ELGIN   G.    FRY        DP*    SRS.  USDA 

C 

DlMLN^alON   X(5u«100)»  A<100) 
REWIND  1 
REWIND  7 


READ 

I  NPUT 

TAPE 

1  < 

(  (X  (  I 

*U)  t 

J=2. 

B  )  •  I 

=  1.4) 

FORMAT ( 12X 

t  7F9. 

=>) 

READ 

I  NPUT 

TAPE 

I  < 

(  (X(  I 

«  J)  . 

J  =  9, 

1  5  )  » 

1  =  1. 

fa) 

READ 

I  NPUT 

TAPE 

1  « 

(  (X(  I 

«  J)  * 

J=16 

*22  )  . 

1  =  1 

.11) 

READ 

I  NPUT 

TAPE 

1  * 

2« 

(  (X(  I 

«  J  )  ♦ 

J  =  23 

♦  29)  « 

1  =  1 

.15) 

READ 

I  NPUT 

TAPE 

1  « 

a* 

(  (X(  I 

»U)  . 

J  =  30 

.35  )  . 

I  =  1 

.18) 

READ 

INPUT 

TAPE 

1  t 

2  • 

{  (X(  I 

.  J)  . 

J  =  37 

*43  }  « 

I  =  1 

.22) 

READ 

I  NPUT 

TAPE 

1  t 

2» 

(  (X  (  I 

.  J)  . 

J  =  44 

«50  )  ♦ 

I  =  1 

♦  25  ) 

READ 

I  NPUT 

TAPE 

1  « 

2< 

{  (X(  I 

«  J  )  « 

J  =  51 

.57  )  . 

I  =  1 

.c:9) 

READ 

I  NPUT 

TAPE 

1  t 

a« 

(  (X(  I 

♦  J)t 

J  =  bb 

.64  )  ♦ 

1  =  1 

.32  ) 

READ 

INPUT 

TAPE 

1  « 

2t 

(  (X(  I 

.J). 

J  =  fcib 

.71  )  . 

I  =  1 

.36) 

REAu 

i  NPUT 

TAPE 

1  « 

2» 

( (X( I « J)» 

J  =  72 

.7b  )  . 

1  =  1 

.39  ) 

READ 

I  NPUT 

TAPE 

1  t 

2  I 

{  (X(  I 

•  J)  • 

J  =  79 

.85)  . 

1  =  1 

.43) 

READ 

INPUT 

TAPE 

1  * 

2» 

(  (X(  I 

»  J)  . 

J  =  ti6 

.92  )  . 

1  =  1 

.46) 

READ 

I  NPUT 

TAPE 

1  . 

2< 

(  (X  (  I 

»  J). 

J  =  93 

.99  )  . 

1  =  1 

.50  ) 

READ 

INPUT 

TAPE 

1  . 

4  t 

(X  I  I  t 

1  00) 

t    1  =  1 

.50  ) 

4   FORMAT ( 12X.F9.5 ) 

DO    14  1=2.100 

DO   6  J=l » 1 
6   A( J )    =  0,0 

K    =  I 

M   =    I    /  2 

DO   a   0=1 .M 

A  (  J  )    =   X  (  J  .  I  ) 

A (K )    =    -X ( J. I ) 

K   =   K   -  1 
8  CONTINUE 

DO    1 O   J=l «  I 
10   A( J)    =   A{ J)    +  4.0 

WRITE   OUTPUT   TAPE   7.    12.    (A(J).  J=l.l) 
12   FORMAT ( 1  OF  10. 5) 
14  CONTINUE 

END  FILE  7 

REWIND  1 

REWIND  7 

Type  i6 
16  format ( 1 ohend  of  job) 

STOP  9999 
END 
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COMPILE        RUN  FORTRAN 
PREPARATION   OF   NORMAL    ORDER    STATISTICS   FOR  DATAPLOTTER 
JOB   NUMBER  920603-0010 
PROGRAM    NUMBER  0911 

WRITTEN   BY   ELGIN   G.    FRY        DP,    SRS,  USDA 

DIMENSION    X(5050).    A(IOO),    P(ll),    CAP(100,10),    SOL ( 1 0 » 1 1 ) ,  N(10) 

REWIND  1 

REWIND  2 

REWIND  3 

REWIND  7 

REWIND  8 

LINE   =  0 

L    =  1 

DO   6  1=2,100 

READ    INPUT    TAPE   7,    2,     (A(J),    J=  1  ,  I  ) 
2    FORMAT ( 1  OF  1 0 .5 ) 

K    =  I*(I-l)/2+l 
NN    =    I    *    (  I    +    1  )    /  2 
J    =  0 

DO   4  M=K,NN 

J    =    J    +  1 
4    X(M)    =    A( J) 
6  CONTINUE 

DO   8  1=2,5050 
8X(I)=X(I)*0.1 

p ( 1  )  =  o»904C  g  feJ-  Cru/ot  r  "  9P99) 

P{2)    =    0.6878  f^^'^ 

P(3)    =  0.6326 

P(4)    =  0.6054 

P<5)    =  0.5751 

P(6)    =  0.5282 

P(7)    =  0.4842 

P(8)    =  0.4000 

P(9)    =  0.3158 

P  (  1  0  )    =  -O.gDOCr 

P  (  1  1  )    =    0. 1674 
10   READ    INPUT    TAPE   8,    12,    STATl,    STAT2,    AREA,    KK  ,  KODE 
12   FORMAT (2A5,F1 0.3, 13, 13 ) 

IF    (KODE-999)  14,50,50 
14    DO    18    1=1 ,KODE 

READ    INPUT    TAPE   8,    16,     (CAP(I,J),  J=l,10) 
16   FORMAT ( 1  OF  1 0.4 ) 
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18  CONTINUE 

READ    INPUT    TAPE   8,    20,    ((SOL(I,J),    J= 1 , 1 1 ) ,  1=1,10) 
20    FORMAT (7X, 1 IFIO.4 ) 

DO   22    I = 1 . 1 0 
22    N (  I  )    =  0 

DO   26    1=1 ,KODE 

DO   26   J= 1  ,  1  0 

IF    (CAP(  I  , J)+9«0  )  24,26,24 
24   N ( J )    =    N { J )    +  1 

26  CONTINUE 

DO   27    1=1 ,KODE 
DO   27   J=l  ,1  0 

27  CAP(I,J)    =    (CAP(I,J)    +   4.0)    *  0.1 
DO   28  1=1,10 

DO   28   J= 1  ,  1  1 

28  SOL{I,J)    =    {SOL(I,J)    +   4.0)    *  0.1 
DO   44    I = 1 , 1 0 

IF    (N( I ) )  44,44,30 
30   K    =   N(I)    *    (N(I)    -    1)    /   2   +  1 
NN=N(I)*(N(I)+l)/2 
J    =  1 

WRITE   OUTPUT    TAPE   L,    32,    STATl,    STAT2 ,    KK ,    I,    X(K),    CAP(J,I),  J 
32   FORMAT (2A5,2 I  3 ,2F8.4 , 1 9x ♦ 1H9, I  8 ) 
K    =   K    +  1 
DO   36  M=K,NN 
J    =    J   +  1 

WRITE   OUTPUT    TAPE   L,    34,    STATl,    STAT2,    KK ,    I,    X(M),    CAP(J,I),  J 
34   FORMAT (2A5,2 I  3 ,2F8. 4 , 1 9x, 1 H8 , IB ) 

36  CONTINUE 
J    =    J   +  1 

WRITE    OUTPUT   TAPE   L,    38,    STATl,    STAT2,    KK ,    I,    X ( NN ) ,    CAP(J-1,I),  J 

38  FORMAT ( 2 A5 , 2  I  3 , 2F8 . 4 , 19x, 1H7, 18 ) 
DO   39   JJ=  1,11 

IF    (SOL(I,JJ))  39,39,41 
41    J    =    J    +  1 
I  I    =  JJ 

WRITE   OUTPUT    TAPE   L,    40,    STATl,    STAT2 ,    KK ,    I,    P ( J J )  ,    SOL (  I  , J J )  ,  J 
40   FORMAT (2A5, 213, 2F8.4,20X, 18 ) 

39  CONTINUE 
J    =    J   +  1 

WRITE   OUTPUT    TAPE   L»    38,    STATl,    STAT2 ,    KK ,    I,    P(II),    SOL (  I  , I  I  ) »  J 
LINE    =   LINE   +   N(I)    +  20 


IF    (LlNE-25000)  44*42* 
42   END   FILE  L 
REWIND  L 
L    =   L   +  1 

LINE  =  0 
44  CONTINUE 

GO  TO  10 
50    REWIND  7 

REWIND  8 

END   FILE  L 

REWIND  L 

TYPE  52 
52   FORMAT ( 1 OHEND   OF  JOB) 

STOP  9999 

END 
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Preparation  of  data  furnished  by 
U.S.  Geological  Survey  (USGS) 


USGS  provides  the  annual  (water  year)  high-flow  volumes  for  periods  of 
1,  3,  T,  15,  30,  60,  90,  120,  150,  l83,  and  2lh  days  for  stream-gage 
records  requested  by  the  Central  Technical  Unit,  SCS . 

The  data  received  from  USGS  are  on  tape  or  ozalid  prints.     If  ozalid 
prints  are  received,  the  IBM  cards  are  key-punched  and  verified  before 
processing.     For  tape,  the  following  programs.  Compile  run  autocoder  by 
P.  Antoinette  Vann  and  format  USGS  data,  FORTRAN  program  No.  Il83,  are 
used  to  punch  cards  directly  from  USGS  tape.     The  header  cards  are  in- 
serted before  listing  for  visual  inspection  and  analysis  on  the  com- 
puter.    Catawba  Creek  near  Catawba,  Va. ,  Station  No.   02-0185.00,  drain- 
age area  =  3^-0  sq.  mi.,  for  19^^-62,  is  a  typical  example. 
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PARAMI9Z1 OlOP 

Q  FORMATS   DATA    FROM  GS 

C  JOB   NUMBER  920603-0030 

(~  PROGRAM   NU^'PER  1183 

n I  WENS  I  ON    X ( 1 O  ) 

1  PEAO    INPUT    TAPE    1,    6«    K,    J*  N 
5   FORMAT ( I?« 18, IX, 18) 

IP    (K-99)  -?,30,30 
30    jp    ( j_g99gq99C>  )  A?*20,20 
32    IF    (K-1  )    34 ,42,42 
34    IF    (N-1  )  36,42,42 
36   PR  I  NT    38 ,    K  »  J 
38   FORMAT ( 1  HI  ,  I  2 ,  I  8 ) 

PUNCH    40,    K,  J 
40   FORMAT (12,18) 

GO   TO  1 
42   BACKSPACE  1 

2  READ  INPUT  TAPE  1,  4,  K,  (X(I)«  1=1' 
4    FORMAT (  I  2 , 7F8. 2 , 3F7, 2 ) 

PRINT    16,    ^  <>     (X(  I  )  ,  1=1*10) 
15   FORMAT ( 1  HO ,  I  2 , 7F8, 1  , 3F7. 1  ) 

pi^lMCH    54, X,     (X(I)»  I=1»10) 
54   F0RMAT( I2,7F8, 1 ,3F7. 1 ) 

GO   TO  1 
20   REWIND  1 
PRINT  22 

22   FORMAT ( 1  HI  ,25X, 1 OHEND   OF   RUN   //// ) 
STOP  999 
END 
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Automatic  data  plottiJag 

The  following  print -out  of  coded  rav  and  computed  data  f or #Catavba  Creek 
and  Smiths  Fork  is  from  tape.    The  automatic  data  plotter^  plots  from  cards 
punched"  from  this  tape  and  a  range  card  made  up  manually. 

1.  Col\imn  1  is  the  station  number  (can  be  any  10  digits  or  symbols). 

2.  Coliimn  2  is  the  set  number  (can  be  10  sets^  only  one  in  this  print-out). 

3.  Column  3  is  the  curve  number  or  field  designation^  which  goes  from  1 
to  10  if  all  fields  are  used.    Both  of  our  examples  contain  10 
coliimns  of  data. 

k.  Column  k  is  the  coded  plotting  position which  is  a  function  of  the 
sample  size  except  as  noted  in  item  6. 

5.  Column  5  gives  the  coded  raw  and  computed  data.    See  item  6. 

6.  Column  6  is  a  code  used  for  line  or  dot  mode  plotting — 9  is  the  code 
on  card  1^  8  on  cards  2  through  19  (sample  size  is  I9  years),  and  7 
a  stop  code  on  card  20.    Cards  21  through  51  are  for  the  11  computed 
points  with  their  respective  coded  plotting  positions  given  in  column 
k  (0.9^85  is  the  0.0  probability  level  and  O.1674  is  the  99. 0  percent 
level) . 

7.  The  last  column  gives  the  card  numbers  for  each  column  or  field.  The 
total  is  the  sample  size  plus  13  (11  probabilities  and  2  stop  cards 
coded  7)«    Catawba  Creek,  N  =  19  (years  of  record)  plus  15,  has  a 
total  of  52  cards  for  each  of  the  10  curves.    Smiths  Fork,  N  =  I6 
(years  of  record)  plus  15,  has  a  total  of  29  cards  for  all  curves 
except  curve  5  for  which  N  =  I5,  a  total  of  28  cards. 

The  ADP  plotting  on  log  normal  probability  paper  for  Catawba  Creek 
follows  the  print -out. 


Dataplotter  Model  5300  manufactured  by  Electronics  Associates,  Inc., 
Long  Branch,  N.  J. 
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FORIRM  program  No.  15^4-6 
"by 

Elgin  G.  Fry 


The  program  is  used  for  plotting  selected  computed  output  of  volume- 
duration  curves  for  dataplotter  on  log-log  paper ^  Volume-Duration- 
Probability  Curves  on  an  Annual  Basis. 

At  the  present  time,  the  selected  computed  values  are  pimched  on  cards 
processed  to  code  the  computed  values  for  plotting  on  the  automatic  data 
plotter.    Modification  of  program  No.  O872  could  give  this  directly  from 
the  main  program,  thus  eliminating  duplication.    The  coding  is  similar 
to  that  used  in  programs  No.  O9IO  and  09II. 
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C 

c 

c 


COMPILE        RUN  FORTRAN 
JOB   NUMBER   920603-9045      PROG  1546 
VOLUME-DURATION   CURVES   FOR  DATAPLOTTER 
WRITTEN   BY   ELGIN    G.    FRY        DP*    SRS,  USDA 


DIMENSION   P(10)»  S0LC4»10) 

REWIND  1 

REWIND  2 

REWIND  3 

REWIND  8 

LINE    =  0 

L    =  1 

P(10)=  0.4438 

P(9)    =  0.4262 

P(8)    =  0.4079 

P(7)    =  0.3954 

P(6)    =    0.3  778 

P(5)    =  0.3477 

P(4)    =  0.3176 

P(3)    =  0.2845 

P(2)    =  0.2477 

P ( 1 )      =    0 .2000 
10    READ    INPUT    TAPE    8,    12*    STAT1»    STAT2*    AREA,    KK »  KODE 
12   FORMAT (2A5 »F1 0.3 »  I  3 ,  I  3 ) 

IF    (KODE-999)  14.50,50 
14   READ    INPUT    TAPE   8,    20,    ((SOL(I,J).    J= 1 , 1 0 ) ,  1=1,4) 
20   FORMAT ( 4X, 1 0F7. 4 ) 

DO   30  1=1,4 

DO    30    J= 1  ,  1  0 

IF    (SOL ( I , J ) )    30 ,30 ,28 
28   SOL (I, J)    =    (LOGXF (SOL ( I , J ) )    +  2.0)    *  0.1 
30  CONTINUE 

J   =  0 

DO    48  1=1,4 
KSW   =  1 
DO    39   JJ= 1 ,10 
IF    (SOL(I, JJ))  39,39,41 
41    J    =    J    +  1 
I  I    =  JJ 

IF    (KSW)  60,60,70 
70   KSW    =  0 

WRITE  OUTPUT  TAPE  L,  72,  STATl,  STAT2,  KK ,  I,  P(JJ),  SOL ( I , J J ) ,  J 
72   FORMAT {2A5,2I 3 ,2F8.4 , 1 9x , 1H9, 18 ) 
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GO    TO  39 

60  WRITE  OUTPUT  TAPE  L«  40,  STATl,  STAT2,  KK ,  I,  P ( J J ) ,  SOL ( I , J J ) i  J 
4  0   FORMAT (2A5 » 2  I  3 , 2F8. 4 , 1 9x , 1 H8 ,  I  8 ) 

39  CONTINUE 
J    =    J    +  1 

WRITE    OUTPUT    TAPE    L,    38,    STATl,    STAT2,    KK ♦     I,    P(II),    SOL (  I  ,  I  I  )  «  J 
38   FORMAT (2A5 ,2  I  3 , 2F8.4 , 1 9x , 1 H7 , IB ) 
48  CONTINUE 

WRITE    OUTPUT    TAPE    L,  80 
80   FORMAT (IH  ) 

LINE    =   LINE    +    J   +  1 

IF    (LINE-25000)  44,42,42 

42  END   FILE  L 
REWIND  L 
KKK    =    L   +  9 

TYPE    43,    LINE,  KKK 

43  FORMAT ( 16 , 16H   RECORDS   ON   UNIT, 13) 
LINE    =  0 

L    =   L    +  1 

44  GO  TO  1 0 
50   REWIND  8 

END   FILE  L 
REWIND   L  . 
KKK    =   L   +  9 
TYPE    43,    LINE,  KKK 
TYPE  52 
52   FORMAT ( 1 GHEND   OF  JOB) 
STOP  9999 
END 
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60    90  120    183  274 
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DURATION    IN  DAYS 


OLUME  -  DURATION  -  PROBABILITY  CURVES  ON  AN  ANNUAL  BASIS 
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C  PERCENTAGE  POINTS  OF  THE   INCOMPLETE  GAMMA-FUNCTION  RATIO  FOR 

C  THE  FOLLOWING  P-VALUESt 

C  P  =  -.95( .05)4 

C  P  =  4.5( .5)38 

C  P  =  39{ 1.0)50 

C  JOB  NO.   920603-3010  PROGRAM  NO.  1319 

C  PROGRAMMED  BY  H.A.RICHARDSON  JULY   13,  1964 

DIMENSION  A(429,5),  B(195,12) 

READ       (   5tl)  ((Ad, J),   J  =   1,   5),   I   =  1,  429) 

1  FORMAT ( lOX, 5E14. 6) 
N  =  0 

DO   2  I   =   1,  39 
DO  2  J  =   1,  11 
K  =   13  -  J 
M  =  5*1 
L  =  M  -  4 
N  =  N  +  1 
IB  =  1 

DO  2   lA  =  L,  M 
B(  IA,K)   =  A(N,IB) 

2  IB  =  IB  +  1 

B(  1,1)  =  -  0.95 
DO   3   I   =  2,  100 

3  B(I , 1)   =  B( I  -   1,1)   +  0.05 
B( 101,1)   =  4.5 

DO  4  I   =   102,  192 

4  B{ 1,1)   =  B( I  -  1,1)   +  0.5 
WRITE  (6,5) 

50F0RMAT(1H1,13X,95HPERCENTAGE  POINTS  OF  THE   INCOMPLETE  GAMMA-FUNCTI 
ION  RATIO  FOR  THE   P-VALUES  INDICATED    IN  COLUMN   1// /I HO , 1 0 X, IHP , 6X, 

2  5H1.000,4X, 5H0.998,4X, 5H0.990,4X,5H0.980,4X ,5H0.960,4X,5H0.900,4X, 

3  5H0.800,4X,5H0.5  00,4X,5H0.200,4X,5H0.100,4X,5HO.OIO} 
WRITE(   9)  ((B(I,J),   J  =   1,    12),    I   =   1,  168) 

WRITE   (6,6)  ((B(I,J),   J  =  1,    12),   I   =  1,  168) 

6  FORMAT(1HO/1HO,4X,12F9.3/(5X,12F9.3) ) 

WRITE(   9)  ((B(I,J),   J  =   1,    12),    I  =  170,    192,  2) 

WRITE   (6,7)  ((B(I,J),   J  =  1,    12),   I   =  170,    192,  2) 

7  F0RMAT(5X,12F9.3) 
WRITE  (6,99999) 

99999  F0RMAT(1H1,20X,10HEND  OF  JOB/lHl) 
STOP 
END 
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c 

JUd    NU    y^UoU:>UUlU  rKUbUof^ 

i' 

t 

c 

HYDRDLuGt    5 1 UuY      bL  b 

j 

L 

JUd    iMUriDCK  V^UOUJ^UUIU 

c 

PROGRAM  NUMBER  0872 

c 

WRITTEN  BY  ELGIN  G.  FRY  DP-SRS-USDA 

1 

c 

INPUT  TAPE   IS  5  AND  2 

c 

OUTPUT   IS  ON  TAPE  3  AND  PRINTER 

1 

c 

WORK  TAPE  ON  4 

1 

c 
c 

END  OF  JOB  CARD — NINES  PUNCHED  IN  COLUMNS 

1-30 

c 
c 

IMPLEMENTED  FOR  SYSTEM  /360  FORTRAN  4  BY 

MERLE  THOMAS  CORP. 

12/65 

ODI MENS  ION  X{ lOOtlO) ,Y( lOOf 10) f lYR(lOO) tXLOG( 100,10) ,CAPX(100,10), 
1IYEAR(100,10),ISTART(10),   lEND(lO),   DAY(IO),   N(10),  OBS(10)t 
2SUMX(10),   SUMXX(IO),  CV(IO),   XM(IO),  GM(IO),   SD(IO),  YX(IO), 
3BETA(10),   BETA2{10)t  GAMMA(IO),   SQGAMdOl,   U(10,ll),   P(  10)  ,CZ(  10), 
4TABLE( 180,12),   XLN(IO),   PK(ll),   XXN(ll),  CHANCE(ll) 
NUM  =  0 

CHANCE(l)   =  0.0 
CHANCE(2)   =  0.2 
CHANCE(3)   =  1.0 
CHANCE(4)   =  2.0 
CHANCE(5)  =  4.0 
CHANCE(6)   =  10.0 
CHANCE(7)   =  20.0 
CHANCE(8)   =  50.0 
CHANCE(9)   =  80.0 
CHANCE(IO)   =  90.0 
CHANCE(ll)   =  99.0 
DAY(l)   =  1.0 
DAY(2)   =  3.0 
DAY(3)   =  7.0 
0AY(4)   =  15.0 
DAY(5)   =  30.0 
DAY(6)   =  60.0 
DAY(7)   =  90.0 
DAY(8)   =  120.0 
DAY(9)   =  183.0 
DAY(IO)   =  274.0 
PK(1)   =  6.015 
PK(2)   =  2.878 
PK(3)   =  2.326 
PK(4)   =  2.054 
PK(5)   =  1.751 
PK(6)   =  1.282 
PK(7)   =  0.842 
PK(8)   =  0.0 
PK(9)   =  -0.842 
PK(IO)   =  -1.282 
PK(ll)   =  -2.326 

C  READ  IN  TABLE  OF  U  VALUES  FOR  VARIOUS  LEVELS  OF  PROBABILITY 

READ(   2)    ( (TABLE( I , J),    J=l,12),  1=1,168) 
READ(   2)    (  (TABLEd  ,J),    J=l,12),    1  =  169,180) 
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C  READ   IN  PARAMETER  CARD 

2  READ   (5,4)    STAT  1 , STAT2 » STAT3 » ARE A , ( I  START ( I ) , I  END ( I ) » 1= 1 , 10 ) 

4  FORMAT(2A4,A2,F10.3t20I3) 
IF    (  ISTART( I)-99 )  5,174,174 

C  READ   IN   DATA  CARDS 

5  K  =    IEND(  1) 

READ   (5,8)    (lYR(I),    ( X (I , J ) , J  =  l  ,  10 )  ,    1  =  1, K) 
8  F0RMAT(I2,7F8.1,3F7.1) 
DO   170  KK=1,10 
Kl   =  ISTART(KK) 
K2  =  lEND(KK) 
K3  =  K2  -  1 
IF    (Kl)  2,2,10 
C  CHANGE   STORAGE  LOCATIONS 

10  DO    12  I=K1,K2 
DO   12  J=l,10 
XLOG(I,J)   =  0.0 
Y( I , J)   =  X(I , J) 

12  IYEAR(I,J)   =  lYR(I) 

C  PLACE  Y  VALUES  ORDERED   HIGH  TO  LOW 

DO   16  J=l,10 

13  DO    16  I=K1,K3 

IF    (Yd, J)   -  Y(I  +  1,J))  14,16,16 

14  HOLD  =  Yd,  J) 

Y(  I  ,  J)    =  Y(  I  +  l,  J  ) 
Y(I+1,J)    =  HOLD 
IHOLD  =    lYEARd  ,  J) 
IYEAR(I,J)    =  IYEAR(I+1,J) 
IYEAR(I+1,J)  =IHOLD 
GO  TO  13 

16  CONTINUE 

C  CALCULATE  LOG  TO  BASE   10  OF  Y  VALUES 

00   18  I=K1,K2 
DO    18  J=l,10 
IF   (Yd, J))  18,18,17 

17  XLOG(I,J)    =  ALOG10(Y( I , J  )) 

18  CONTINUE 

C  CALCULATE  CAP(X)  VALUES 

IF    (AREA)  19,19,22 

19  DO   20  I=K1,K2 
DO   20  J=l,10 

20  CAPX(  I  ,  J)   =  Yd,  J) 
GO  TO  26 

22  CC  =  0.03719008  /  AREA 

DO   24  I=K1,K2 

DO   24  J=l,10 
24  CAPX(I,J)    =  Yd, J)    *  DAY(J)    *  CC 
C  WRITE  HEADER  LINES 

26  WRITE( 6,28)STAT1,STAT2,STAT3,AREA,KK 

280F0RMAT( IHl, 10X,61HHIGHEST  MEAN  DISCHARGE  GAMMA  FUNCTION  PARAMETERS 
1   FOR  STATION   , 2A4, A2,4X, 6HAREA  = , Fl 0 . 3 , 5X, 3HSET , I  3) 
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WRITE(6,30) 

300F0RMAT(1H0,115HYR  Xd)        LOG  X(l)       CAP  X(ll  YR  X( 

12)       LOG  X(2)       CAP  X(2)  YR  X(3)       LOG  X(3)       CAP  X(3)  ) 

WRITE  DATA  VALUES 
DO  32  I=K1,K2 

320WRITE   (6,34)    I YE AR { I , 1 ) , Y( I , 1 ) , XLOG ( 1 , 1 ) t CAPX ( I , 1 )  , 

1IYEAR(I,2),   Yd, 2),   XL0G(I,2),   CAPX(I,2),   IYEAR(I,3),  Yd, 3), 
2XL0G(I,3),  CAPX(I,3) 

34  FORMAT   dX ,  I  2,  F 1 2  .4,  F 1 0 .  6,  F  10 . 4 ,  2  d  8  ,  Fl 2  .  4,  F 10.  6 ,  FIO  .4)  ) 

WRITE(6,36) 

360F0RMAT(1H0,115HYR  X(4)        LOG  X(4)       CAP  X(4)  YR  X( 

15)       LOG  X(5)       CAP  X(5)  YR  X(6)       LOG  X(6)       CAP  X{6)  ) 

WRITE  DATA  VALUES 
DO  38  I=K1,K2 

38CWRITE(6,34)    IYEAR(I,4),   Y( I , 4) , XLOG ( I ,4 ) ,  CAPX{I,4), 

1IYEAR(I,5),   Yd, 5),   XL0G(I,5),   CAPX(I,5),   IYEAR(I,6),   Yd, 6), 
2XL0G(I,6),  CAPX(I,6) 
WRITE(6,40) 

40CFORMAT(1HO,115HYR  X(7)        LOG  X(7)       CAP  X(7)  YR  X( 

18)        LOG  X(8)       CAP  X(8)  YR  X(9)       LOG  X(9)       CAP  X(9)  ) 

WRITE  DATA  VALUES 
DO  42  I=K1,K2 

42CWRITE(6,34)    IYEAR(I,7),   Yd, 7),   XL0Gd,7),  CAPXd,7), 

lIYEARd,8),   Yd, 8),   XL0G(I,8),   CAPXd.B),   IYEAR(I,9>,  Yd, 9), 
2XL0G(I,9),  CAPX(I,9) 
WRITE{6,44) 

44  FORMAT  (1H0,35HYR  XdO)     LOG  XdO)   CAP  X(IO)  ) 

DO  46  I=K1,K2 

46  WRITE(6,34)    IYEAR(  I  ,  10  )  ,  Yd  ,  10)  ,   XLOG(  I  ,  10  )  ,  CAPX  ( I  ,  10  ) 

DO  41  I=K1,K2 

DO  41  J=l,10 

IF   (CAPX(I,J))  55,55,43 
55  CAPXd,J)   =  -9.0 

GO  TO  41 
43  CAPXd,J)   =AL0G10(CAPX  ( I  ,J)  ) 
41  CONTINUE 

KODE  =  K2  -  Kl  +  1 

WRITE( 3)  STAT1,STAT2,STAT3,AREA,KK,K0DE 

WRITE   (3)    ({CAPX(I,J),   J=l,10),  I=K1,K2) 

CALCULATE  VALUE  OF  N  FOR  EACH  DAY 

DO  48  1=1,10 
48  N( I )  =  0 

DO   52  I=K1,K2 

DO  52  J=l,10 

IF   (Yd, J))   52,52,  50 
50  N( J)   =  N(J)   +  1 
52  CONTINUE 

WRITE(6,53)AREA  c.  .-^tc 

530FORMAT(1H1,30X,66HVOLUME-DURATION-PROBABIL1TY  ANALYSIS  FOR  SELECTE 
10  WATERSHEDS  /1H0,14HNAME  OF  STREAM ,44X , 13HGAGE  LOCATION/ 

21H0,8HUSGS  NO. , 1 5X,7HCTU  NO. , 1 1 X , 16HPERI0D  OF  RECORD, 32X, 
315HDRAINAGE  AREA  =,F8.3,6H  SQ  MI   /IHO , 20HDURAT I  ON  IN  DAYS  1,10X, 
41H3,10X,1H7,9X,2H15,9X,2H30,9X,2H60,9X,2H90,8X,3H120,8X,3H183,8X, 

53H274  ) 
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DO  54  1=1,10 
54  OBS(  I  )  =  N(  I  ) 

WRITE(6,56)    (N(I),  1=1,10) 
56   FORMATdOHON  ,10111) 
C  CALCULATE  MEAN,  CV 

DO  5  8  1=1,10 
XM( I  )   =  0.0 

cv; I )  =  0.0 

CZ(  I  )   =  0.0 
SD( I  )  =  0.0 
XLN(  I )   =  0.0 
GM( I  )  =  0.0 
SUMXd  )  =  0.0 
58   SUMXX(  I  )   =  0.0 
DO  60  I=K1,K2 
DO  60  J=l,10 

SUMX( J)    =   SUMX( J)   +  Y( I , J) 

60  SUMXX(J)    =   SUMXX(J)   +  Yd, J)    *  Yd, J) 
DO  62  J=l,10 

IF    (OBSCJ))  62,62,61 

61  XM(J)   =   SUMX(J)    /  OBS{J) 

OCV(J)  =  ((OBS(J)  *  OBS{J)  *  SUMXX(J)  -  OBS(J)  *  SUMX(J)  *  SUMX(J)) 
1/   ((OBS(J)   -  1.0)   *  SUMX(J)    *  SUMX(J)))   **  0.5 

62  CONTINUE 

C  WRITE  MEANS,  CV,   AND  SD 

WRITE(6,64)    (CVd),   1  =  1,10) 
64   FORMATdOH  CV  ,10F11.5) 

WRITE(6,66)    (XM(I),  1=1,10) 
66  FORMATdOH   X  BAR  ,10F11.4) 
C  CONVERT  LOG  AREA  TO  LOG  OF   X  TO  BASE  E 

00  70  I=K1,K2 
DO  70  J=l,10 
IF    (Yd, J))  70,70,69 

69  XLOG(I,J)   =  ALOG(Y(I,J)) 

70  CONTINUE 

C  CALCULATE   GEOMETRIC  MEAN  AND  C(Z) 

DO  72  1=1,10 
SUMXX( I  )  =  0.0 
72   SUMXd  )   =  0.0 
DO   74  I=K1,K2 
DO  74  J=l,10 

SUMXX(J)   =  SUMXX(J)   +  XLOG(I,J)   *  XLOG(I,J) 
74  SUMX{J)   =  SUMX(J)   +  XLOG(I,J) 

DO  76  J=l,10 

IF   (OBS(J))  76,76,75 
750CZ(J)   =   ((OBS(J)   *  OBS{J)   ♦  SUMXX(J)   -OBS(J)   ♦  SUMX(J)  *SUMX(J)) 
1/   ((OBS(J)   -   1.0)   *  SUMX(J)    *  SUMX(J)))   **  0.5 

SD(J)   =  SUMX(J)    /  OBS(J) 

XLN{J)   =  SD(J) 

CZ( J)   =  CZ( J)   *  ABS   (XLN(J) ) 
GM(J)  =  EXP  (SD(J)) 
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76  CONTINUE 
C  WRITE  LOG  GM  AND  THE  GM 

WRITE(6,78)    (SD(I),  1=1,10) 
78  FORMATdOH  LN  G  ,10F11.6) 

WRITE(6,80)    (GM(I),  1=1,10) 

80  FORMATdOH  G  ,10F11.4) 

C  CALCULATE  AND  WRITE  RATIO  OF  X  BAR  /  G 

DO   82  1=1,10 
SD( I  )  =  0.0 
IF    (OBS(I))  82,82,81 

81  SD( I  )  =  XM( n   /  GM( I ) 

82  CONTINUE 

WRITE(6,84)  (SD(I),I=1,10) 

84  FORMATdOH  RATIO  X/G,10F11.5) 
C           CALCULATE  AND  WRITE  CV*  AND  Y 

DO   86  1=1,10 

YX( I  )   =  0.0 

SUMXd  )   =  0.0 

IF    (OBS(I))  86,86,85 

85  SUMXd)   =   (SD(I)   ♦  SD(I)  -   1.0)   **  0.5 
YX( I )   =  ALOG(SD( I ) ) 

86  CONTINUE 

WRITE(6,88)    (SUMXd),    1  =  1,10) 
88  FORMATdOH  CV*  ,10F11.5) 

WRITE(6,90)    (YX(I),  1=1,10) 
90  FORMATdOH  Y  ,10F11.6) 
C  CALCULATE  LN  S.D.   AND  LN  S.D.*    (USE   BETA  AREA) 

DO  240  1=1,10 

BETA(  I)   =  0.0 

IF    (OBS(I))  240,240,241 

241  BETA(I)   =   (YX(I)   *  2.0)   **  0.5 
240  CONTINUE 

WRITE(6,242)    ( BETA d ) , 1  =  1 , 10 ) 

242  FORMATdOH  LN  S.D.  ,10F11.6) 
DO  244   1  =  1,10 

IF    (OBS(I))  244,244,243 

243  BETAd)   =  BETA(I  )   *   (OBSd)    /   (OBS(I)   -  1.0))   **  0.5 

244  CONTINUE 

WRITE(6, 246)  (BETAd  ),    1  =  1,10) 

246  FORMATdOH  LN  S.D.*  ,10F11.6) 
WRITE(  6,245)    (CZd),    1  =  1,10) 

245  FORMATdOH  LN  S.D.**,  lOF  11.  6) 

C  CALCULATE  BETA  AND  GAMMA  AND  SQR  ROOT  OF  GAMMA 

DO  99  1=1,10 
IF    (OBS(I))  247,247,248 

247  GAMMA( I )   =  0.0 
BETA( I )   =  0.0 
SQGAM( I )   =  0.0 
GO  TO  99 

248  IF   ( YX( I )-0  .5772)  92,94,94 

920GAMMA(I)   =   (1.0  /  YXd))   *   (0.5000876  +  0.1648852  *  YX(I) 
1-  0.0544274  *  YXd)   *  YXd)) 
GO  TO  98 
94   IF    (YXd)   -   17.0)  96,96,170 
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960GAMMA(I)   =    (8.898919  +   9.059950  *  YX(I)   +  0.9775373   *  YX(I)  * 

lYX(I))    /   (YX(I)    *   (17.79728  +   11.968477  *  YX(I)    +  YX(I)    *  YX(I))) 

98  BETAd)   =  XM(I)    /  GAMMA(I) 
SQGAM(I)    =   GAMMA(I)   **  0.5 

99  CONTINUE 
GO  TO  104 

C  WRITE  GAMMA,    SQR   ROOT  OF  GAMMA   AND  BETA 

104  WRITE(6,106)    (GAMMA(I),  1=1,10) 
106   FORMATdOH  GAMMA  ,10F11.5) 

WRITE{6,108)    (SQGAM(I),  1=1,10) 
108   FORMATdOH   SQR  G  AMMA  ,  1  OF  1 1 .  5  ) 

WRITE(6,11C)    (BETAd),  1=1,10) 
110   FORMATdOH   BETA  ,10F11.5) 
C  CALCULATE  AND  WRITE  BETA  STAR  (PRIME) 

DO   116  1=1,10 

IF    (AREA)  114,114,112 

112  IF    (OBS(I))  114,114,113 

113  BETA2(I)   =  DAY(I)   *  0.03719  *  BETAd)   /  AREA 
GO  TO  116 

114  BETA2( I )   =  BETA(  I  ) 
116  CONTINUE 

WRITE(6,118)    (BETA2(I ) ,1=1,10) 
118  FORMATdOH   BETA*  ,10F11.5) 
C  SUMX  =  BETA  X  SQR  ROOT  GAMMA 

C  SD  =  BETA*  X  SQR  ROOT  GAMMA 

DO   120  1=1,10 

SUMXd)    =   BETAd)    *  SQGAM(I) 
120  SD(I)  =  BETA2(I)   *  SQGAMd) 

WRITE(6,122)    (SUMXd),  1=1,10) 
122  FORMATdOH  B  *  S  (  GA  )  ,  1  OF  11 . 4) 

WRITE{6,124)    (SDd),  1=1,10) 
124  FORMATdOH  B*    (  S  (  GA  )  ,  lOF  11 .  6 ) 

DO    126  1=1,10 

126  P ( I )   =   GAMMA( I  )    -  1.0 

C  CALCULATE  AND  STORE  U  ♦  B  *  S(GA) 

DO  127  1=1,10 
DO   127  J=l,ll 

127  U( I  ,  J)   =  0.0 
DO   152  J=l,10 

IF   (OBS(J))  152,152,129 

129  IF   (P(J)-50.0)  128,128,152 

128  DO   130  1=2,180 

IF    (P( J)-TABLE(I,1) )  131,144,130 

130  CONTINUE 
GO  TO  152 

131  IF    ( (P(J)-TABLE( 1-1,1) )   - ( TABLE ( I , 1 ) -P ( J )  )  )  132,132,133 

132  PROB  =   (P(J)  -  TABLE( 1-1,1) )   /   (TABLE(I,1)   -  TABLE( I-l , 1 )  ) 
SW  =  0.0 

GO  TO  134 

133  PROB  =   (TABLE( I , 1)-P( J ) )   /   (TABLE(I,1)   -  TABL E ( I-l , 1 ) ) 
SW  =  1.0 
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134  00   142  M=l,ll 

IF   (TABLE( I-1,M+1)  -  T ABLE ( I , M+1 ) )  138,136,136 

136  IF   (SW)  137,135,137 

137  PROBT=  TABLE(I,M+1)   +   ( ( TABL E ( I -1 ,M+1 )   -  T ABL E ( I , M+ 1 ) )   *  PROB) 
GO  TO  140 

135  PROBT  =  TABLE( I-1,M+1)   -  ( ( TABL E ( I- 1 , M+ 1 )-TABLE (I ,M+1) )   *  PROB) 
GO  TO  140 

138  IF   (SW)  139,139,141 

139  PROBT=  TABLE( I-1,M+1)   +   ( (TABLE ( I ,M+1)  -  TABL E ( I-l , M+ 1 ) )  *  PROB) 
GO  TO  140 

141  PROBT  =  TABLE(I,M+1)  -   { (T ABLE( I , M+1 )   -  TABLE ( I-l, M+1 ) )   ♦  PROB) 

140  U(J,M)   =  PROBT*  SD(J) 

142  CONTINUE 
GO  TO  152 

144  DO   146  M=l,ll 

146  U(J,M)   =  TABLE(I,M+1)   *  SD(J) 

152  CONTINUE 

WRITE  U  *  B  *  S(GA)  VALUES 
WRITE(6,148) 

148  FORMAT( IHO, 14HPERCENT  CH ANCE ,43X , 9HU ( B*S ( GA ) ) 

WRITE(6,154)    (CHANCEd),    (U(J,I),   J  =  l,10),   1  =  1,11) 

154  FORMATdH   ,  F9  .  1 ,  lOF  1 1 .  4 ) 
WRITE( 6, 161  ) 

1610FORMAT(120H0N  =  STATION  YEARS,  CV  =  COEFFICIENT  OF  VARIATION  BY  ME 
ITHOD  OF  MOMENTS,  X  BAR  =  ARITHMETIC  MEAN  IN  SECOND-FEET  PER  DAY,  / 
2116H  LN  G  =  NATURAL  LOGARITHUM  OF  THE  GEOMETRIC  MEAN(MEAN  OF  THE  N 
3ATURAL  LOGARITHMS  OF  THE  ORIGINAL  DATA),  G  =  GEOMETRIC/  120H  MEAN 
4IN  SECOND-FEET  PER  DAY,  RATIO  X/G  =  RATIO  OF  ARITHMETIC  MEAN  TO  GE 
50METRIC  MEAN,  CV*  =  COEFFICIENT  OF  VARIATION  BY  /120H  METHOD  OF  MA 
6XIMUM  LIKELIHOOD,  Y  =  NATURAL  LOGARITHM  OF  THE  RATIO  OF  THE  ARITHM 
7ETIC  MEAN  TO  THE  GEOMETRIC  MEAN,  LN  S.D.  /114H  AND  LN  S.D.*  =  STAN 
8DARD  DEVIATION  OF  THE  NATURAL  LOGARITHMS  FOR  MAXIMUM  LIKELIHOOD,  B 
9IASE0  AND  UNBIASED  RESPECT-  ) 
WRITE(6, 185) 

1850F0RMAT( 120H  IVELY,  GAMMA  =  SHAPE  STATISTIC,  SQR  GAMMA  =  SQUARE  ROO 
IT  OF  GAMMA,  BETA  AND  BETA*  =  RATIO  OF  ARITHMETIC  MEAN  TO  GAMMA  IN/ 
2116H  SECOND-FEET  PER  DAY  AND  VOLUME  IN  INCHES  RESPECTIVELY,  B  *  S( 
3GA)    AND  B*   (S(GA)   =  BETA  TIMES  THE  SQUARE  ROOT  OF  /120H  GAMMA 

4(STANDARD  DEVIATION)  IN  SECOND-FEET  PER  DAY  AND  VOLUME  IN  INCHES  R 
5ESPECTIVELY,  U(3*S(GA)  =  VOLUME  IN  INCHES  FOR  /  51H  SELECTED  DUR 
6ATI0NS  AND  PERCENT  CHANCE  AS   INDICATED  /  IHl  ) 

DO   155  1=1,10 

DO   155  J=l,ll 

IF   (U(I,J))  151,151,153 
151  U{ I , J)   =  -9.0 

GO  TO  155 

153  U(I,J)   =ALOG10(U( I,J) ) 

155  CONTINUE 

DO   166  1=1,10 
IF    (OBS(I))  156,156,157 
157   IF    (GAMMA( I )-51.0)  156,156,158 

156  WRITE(3)  (U(I,J),  J=l,ll) 
GO  TO  166 

CALCULATE  LOG-NORMAL  SOLUTION 


158  DO   160  J=l,ll 

XXN(J)    =   EXP    (XLN(I)   +   PK(J)    *  ABS  (CZ(I))) 

IF    (AREA)  160,160,183 
183   XXN(J)    =   XXN(J)   *  DAY(I)    *  0.03719   /  AREA 
160  CONTINUE 

WRITE(6,162)    I,    (XXN(J),  J=l,ll) 

162  FORMAK IHO, I2,4X,11F10.6) 
00    163  J=l,ll 

IF    (XXN(J))  171,171,165 

171  XXN(J)    =  -9.0 
GO   TO  163 

165  XXN(J)    =AL0G10(XXN( J)  ) 

163  CONTINUE 

WRITE(3)  (XXN(J),  J=l,ll) 

166  CONTINUE 

170  NUM  =   NUM  +    (4  *  K2  +  37) 
IF   (NUM-20000)  2,172,172 

172  NUM=0 
GO  TO  2 

174  CONTINUE 

KODE  =  999 

WRITE   (3)  STAT1,STAT2,STAT3, AREA, KK, KODE 

WRITE  (6,99999) 
99999  FORMAT( 1H1,20X,10HEND   OF  JOB/lHl) 
STOP  9999 
END 
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C  CONVERSION  OF  NORMAL  ORDER  STATISTICS  ' 

C  JOB  NUMBER  920603-0010  ' 

C  PROGRAM  NUMBER  0910 

C  WRITTEN  BY  ELGIN  G.   FRY       DP,   SRS,  USDA 

C  IMPLEMENTED  FOR  SYSTEM  /360  FORTRAN  4  BY  MERLE  THOMAS  CORP.  12/65 
C 

DIMENSION  X(50,100),  A(ICO) 
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READ 

(It 

2) 

(  (X(  I 

tJ)t 

J= 

58,64), 

1  =  1 

,32) 

READ 

(1» 

2) 

(  (X(  I 

t  J)t 

J= 

65,71) , 

1=1 

,36) 

READ 

(It 

2) 

(  (X(  I 

t  J)t 

J= 

72,78), 

1  =  1 

,39) 

READ 

(It 

2) 

(  (X(  I 

t  J)t 

J= 

79,85), 

1  =  1 

,43) 

READ 

(It 

2) 

(  (X(  I 

t  J)t 

J= 

86,92) , 

1=1 

,46) 

READ 

(  It 

2) 

(  (X(  I 

t  J)t 

J= 

93,99) , 

1=1 

,50) 

READ 

(It 

4) 

(X(I, 

100), 

I 

=1,50) 

DO   14  1=2,100 
DO  6  J=1,I 
6  A(J)  =  0.0 
K  =  I 
M  =   I   /  2 
DO   8  J=1,M 
A(J)   =  X(J,I) 
A(K)  =  -X( J, I ) 
K  =  K  -  1 
8  CONTINUE 
DO   10  J=1,I 
10  A(J)   =  A(J)   +  4.0 

WRITE   (4)  (A(J)    ,  J=1,I) 

14  CONTINUE 

WRITE  (6,99999) 
99999  F0RMAT(1H1,20X,10HEND  OF  JOB/lHl) 
STOP  9999 
4  F0RMAT(12X,F9.5) 
2  F0RMAT(12X,7F9.5) 
END 


! 
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r  i?^!;^"^'^^^^  SYSTEM  /360  FORTRAN  4     BY  MERLE  THOMAS  CORP 

C  PREPARATION  OF  NORMAL   ORDER  STATISTICS  FOR  DATAPLOTTER 

C  JOB  NUMBER  920603-0010  UAiAPLUIItR 

C  PROGRAM  NUMBER  0911 

C  WRITTEN  BY  ELGIN  G.   FRY       DP,   SRS,  USDA 

DoTllT^lOo'''''^'  CAP(100,10),   SOL(10,11),  N(IO) 

READ  (2)        {A(J),  J=1,I) 

K=I*(I-l)/2+l 

NN  =  I   *   (I   +   1)   /  2 

J  =  0 

DO  4  M=K,NN 
J  =   J   +  1 
4  X(M)    =  A(J) 
6  CONTINUE 

DO  8  1=2,5050 
8   X( I )   =  X( I )   *  0.1 

P(l)   =  0.9999 

P(2)   =  0.6878 

P(3)   =  0.6326 

P(4)    =  0.6054 

P(5)   =  0.5751 

P{6)   =  0.5282 

P(7)   =  0.4842 

P(8)   =  0.4000 

P(9)   =  0.3158 

P(10)  =0.2718 

P(  11)   =  0.1674 

^I'^^.lll         ^^"^^^  STAT3,   AREA,   KK.  KODE 

IF   (KODE-999)  14,50,50 

14  READ   (3)        ((CAP{I,J),    J=l,10),    1=1, KODE) 
DO   15  1=1,10 

15  READ   (3)  (SOL(I,J)  ,  J=iai) 
DO  22  1=1,10 

22  N(  I  )  =  0 

DO  26   1=1, KODE 
DO  26  J=l,10 

IF   (CAP( I,J)+9.0)  24,27,24 
24  N(J)   =  N(J)   +  1 

27  CAP{I,J)   =   (CAP(I,J)   +  4.0)   *  0.1 
26  CONTINUE 

DO  28  1=1,10 
DO  28  J=l,ll 

28  SOL(I,J)   =   (SOL(I,J)   +  4.0)   *  0.1 
DO  44  1=1,10 

IF    (N(I))  44,44,30 
30  K  =  N(  I)   *   (Nd)   -  1)    /  2  +  1 
NN  =  N(I)   ♦   (N(I)   +  1)    /  2 
J  =  1 

WRITE   (4,90)   STATl,   STAT2 , STAT3 ,KK ,    I,   X(K),   CAP(J,I),  j 
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90  FORMAT   ( IHV, 2 A4 , A2 , 2 1  3 , 2 F8. 4, 1 9X , 1H9, I  8 ) 
K  =  K  +  1 

DO  36  M=K,NN 
J  =  J  +  1 

WRITE   (4,91)   STATl,   ST  AT  2, STAT3 , KK, I » X ( M ) , C AP ( J , I ) , J 

91  FORMAT   ( lHV,2A4t A2,2I3,2F8,4, 19X,1H8,I8) 

36  CONTINUE 
J  =  J  +  1 

WRITE   (4,92)   STATlt   STAT2tSTAT3,   KK,    I,   X(NN),   CAP(J-1,I),  J 

92  FORMAT   ( IHV » 2 A4, A2 » 21 3 , 2F8. 4, 19X , 1H7, I  8 ) 
DO  39  JJ=1,11 

IF   (SOL(I,JJ))  39,39,41 
41   J  =  J  +  1 
II   =  JJ 

WRITE   (4,93)   STATl,   ST AT2 , STAT3 , KK , I , P ( J J ) , SOL ( I , JJ ) , J 

93  FORMAT   { IH V , 2 A4 , A2 , 2 1 3 , 2F8. 4, 20X , I  8 ) 

39  CONTINUE 
J  =  J  +  1 

WRITE   (4,92)   STATl,   STAT2 , STAT3 , KK, I , P( 1 1 ) , SOL ( I , 1 1 ) , J 

44  CONTINUE 
GO  TO  10 
50  CONTINUE 

WRITE  (6,99999) 
99999   FORMAT( 1H1,20X,10HEND   OF  JOB/lHl) 
STOP  9999 
END 


000010   IDENTIFICATION  DIVISION, 
000020   PROGRAM-ID.  'AOIl^l'. 

000030  AUTHOR.  BRANDON  APPLIED  SYSTEMS,    INC,   G  POPKIN. 

00004C   REMARKS.   NOTE-IT — FDA114110  CANNOT   BE  BLOCKED. 
000050   ENVIRONMENT  DIVISION. 
00006C  CONFIGURATION  SECTION. 
000070  SOURCE-COMPUTER.    IBM-360  G40. 
000080  OBJECT-COMPUTER.   IBM-36C  GAO. 
000C90  INPUT-OUTPUT  SECTION. 
000100  FILE-CONTROL. 

000110  SELECT  MESSAGE  ASSIGN  TO  'MSGOUT*  UTILITY. 

000120  SELECT  A114151   ASSIGN  TO   'OUTPUTl'  UTILITY. 

000130  SELECT  A114110  ASSIGN  TO   'INPUTI'  UTILITY. 

000140  DATA  DIVISION. 
000150  FILE  SECTION. 

000160  FD     MESSAGE  RECORDING   F  LABEL  RECORD  STANDARD  DATA  RECORD 
000170  01     MES-TM   PICTURE  X(160). 
000180  FD  A114110 

000190  BLOCK  CONTAINS  0001  RECORDS 

000200  RECORD  CONTAINS  0165  CHARACTERS 

001000  LABEL   RECORDS   ARE  STANDARD 

OOlOlC  DATA  RECORDS  ARE  INPUTX 

00102C  RECORDING  MODE    IS  F. 

001030  01  INPUTX. 

00104C  02  INPUT-XXXXXGENX. 

001050  08     STATNO  PICTURE  S9(011). 

00106C  08     FILLER  PICTURE  X(Oll). 

00107C  08     WATERYR  PICTURE  S9(011). 

001080  08     HIDAY  PICTURE  S9(011). 

001090  08     HIDAY3   PICTURE  S9(011). 

001100  08     HIDAY7  PICTURE  S9(011). 

OOlllC  08     HIDAY15   PICTURE  S9(011). 

001120  08     HI0AY30   PICTURE  S9(011). 

001130  08     HIDAY60  PICTURE  S9(011). 

001140  08     HIDAY90  PICTURE  S9(011). 

00115C  08     HIDAY120  PICTURE  S9(011). 

001160  08     FILLER  PICTURE  X(Oll). 

001170  08     HIDAY183  PICTURE  S9(011). 

00118C  08     HIDAY274  PICTURE  S9(011). 

00119C  08     1GENDD0012  REDEFINES  HIDAY274. 

001200  10     FILLER  PICTURE  X(004). 

002000  10     HIDAY274-XXXXX0083   PICTURE  S9{003). 

002010  10     HIDAY274-XXXXX0141   PICTURE  S9(004). 

002020  02     FILLER  PICTURE  X(ll). 

002030   FD  A114151 

002040  BLOCK  CONTAINS  0001  RECORDS 

002050  RECORD  CONTAINS  0080  CHARACTERS 

002060  LABEL   RECORDS   ARE  STANDARD 

002070  DATA  RECORDS  ARE  OUTPUTX 


MES-TM. 
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Wc.  uou 

RECORDING  MODE 

IS  F. 

1 1 A 1 

KJU  c.  U  7  u 

01 

OUTPUTX. 

02 

OUTPUT-XXXXGENX. 

1  1  1 

1141 

u  u  ^  i  1  u 

08 

FILLER  PICTURE  X(080) 

• 

1  1  A  1 
1  iHl 

o  n    1  9  n 

WORKING 

-STORAGE  SECTION. 

1141 

UUc,  ldKJ 

n  1 

U  1 

IIWX. 

1  1  A  1 

1 141 

n  no  1  An 

02 

2IW99  PICTURE  9(10) 

• 

1  1  A  1 

1141 

r\r\o  ^  «^ n 
UU^  i  OV 

02 

FILLER  REDEFINES  2IW99. 

1  1  A  1 
1  IHI 

yjyjc.  L  ou 

03 

FILLER  PICTURE  99. 

1  1  A  1 
1  1  Hi 

C\f\0  1  ICS 

03 

1IW99  PICTURE   9(4) . 

1  1  A  1 
1  IHI 

f\c\o  ^  DO 
uu  ^ 1 ou 

03 

0IW99  PICTURE  9(4). 

1  1  A 1 

r\     o  1  on 

n  1 

U  X 

lELSW. 

1  1  A  1 
1  iHl 

nno  o  n A 

n? 

\J  c. 

lGENERELSW-001 

PICTURE 

y    V  ALUt    U . 

1  1  A 1 
1  l'»l 

UU  :5  UUU 

n  7 

1  Fl  <;un 7 

let-  OwtKiC 

PT  rXIIRF 

q 

VAI  IIF 

U  . 

1  1  A 1 
1  It-i 

UU  J  U  i  u 

u  c. 

J.  C  l_  O  W  v  ^ 

P  I  TTI IR  F 

q 

\l L\  IIF 

u. 

1  1  A  1 
1  If-l 

0? 

1  Fl  <s W04 

PTCTURF 

9 

VAI  IIF 

V  «  L  "L»  L 

n 

1  1  Al 
1  IHl 

U  U J  U3  u 

X  c  i_  o  n  w  ^ 

9 

M  L\  IIF 

V  M 1-  C 

n 

U  . 

X  1*H 

n? 

1  Fl  <;unA 

P  T  r  Tl  IR  F 

9 

M  L\  IIF 

Ua 

1  1  Al 
X  X**l 

yjyj  D  UP  u 

OP 

1  Fl  <;wn7 

J.  1^  1-  o  n  v  1 

P  T  r  Tl  IR  F 

q 

V&l  1  IF 

n 

1  1  A.1 

n H"^  HAH 

0? 

1  Fl  ^UDR 

P  I  PT  1  IR  F 

q 

\/ Al  GIF 

V  ML  UC 

n 
u  . 

1  l*f  X 

U  U  3  U  1  u 

1  Fl  <\WnQ 

P  T  TTIIRF 

9 

V  Al  IIF 

V  M  l_  C 

u  • 

1  1  A  1 
X  X*tX 

uu  ^  uo  u 

1  F 1  <\  U 1  n 
X  c  L  o  n  X  \J 

P T  r THR  F 

q 

U  Al  IIF 

V  H  L  W  ^ 

n 
u  . 

X  X  T-X 

U  U  J  u 

1  F  1  ^  U  1  1 
X  C  L.  O  n  X  X 

P  T  TTIIR  F 

q 

VAI  IIF 

V  ML.  U C 

u  . 

1  1  Al 
X  X*tX 

UUO 1 UU 

1  Fl  "nUI  7 
X  C  l_  On  X  c 

p  T  r  Tl  IR  F 

q 

\/ Al  IIF 

V  ML  UC 

u  • 

1  1  Al 
X  XHX 

n  n*^  1  1  n 
u  u  ^  i  i  u 

1  Fl 

X  1^  L.  O  T1  X  ^ 

PT  PTURF 

9 

VALUE 

n 

1  1  Al 
X  X*rX 

nn    1  on 
UU  3  i  ^u 

07 

1  Fl  <IU1  U 
X  c  l_  o  W  X  " 

P  T  r  Tl  IR  F 

q 

U Al  IIF 

V  M  L  UC 

n 
u  . 

1  1  Al 
X  X  fx 

n  n  "2  1  "an 
UU -5  i ^ U 

n? 

u  ^ 

1ELSW15 

PICTURE 

9 

VAI  ItF 
y  w  i_  u  t 

u  • 

1  1  Al 
X  X  *f  X 

nn'j  1  />  n 
UUO  1  *f  u 

1ELSW16 

PICTURE 

9 

VALUE 

n 

1  1  Al 
X  X  T-X 

n  n'2  1  c;  n 
U  U  3  J.!3  U 

fl  7 

1ELSW17 

PICTURE 

9 

VALUE 

u  . 

1  1  Al 
X  X  H  X 

n  n'^  1  An 
u  U^ 1 OU 

02 

1ELSW18 

PICTURE 

9 

VALUE 

u  • 

1  1  Al 
X  X  *T  X 

n  nra  1  7n 
U  U^  i  f  U 

02 

1ELSW19 

PICTURE 

9 

VALUE 

u  . 

1  1  A 1 
X  X  *t  X 

nn'a  1  on 
uu^  i ou 

02 

1ELSW20 

PICTURE 

9 

VALUE 

u  • 

1  1  A 1 
X  If  X 

n  n  "3  1  on 

02 

1ELSW21 

PICTURE 

9 

VALUE 

U  . 

1  1  A  1 
X  X  H  X 

on 3  on o 

UU^  C\)\J 

02 

1ELSW22 

PICTURE 

9 

VALUE 

u  . 

1  1  Al 
X  X  *tX 

UU*f  UUU 

02 

1ELSW23 

PICTURE 

9 

VALUE 

n 

1  1  Al 
X  X  *tX 

n  nZt  n  1  n 
UUH  U  i  u 

02 

1ELSW24 

PICTURE 

9 

VALUE 

u  . 

1  1  Al 
X  X*tX 

n  nA  n  o  n 
U  U*f  Uii  u 

02 

1ELSW25 

PICTURE 

9 

VALUE 

n 

u  . 

1  1  Al 
X  X*tX 

no  A  na  n 

02 

1ELSW26 

PICTURE 

9 

VALUE 

n 

1141 

X  X  ^  X 

o  n A  nAn 
U  U*fr  U*f  u 

02 

1ELSW27 

PICTURE 

9 

VALUE 

u  • 

1  1  Al 
X  X  X 

o n A  OR  n 

U  U*T  U-3  u 

02 

1ELSW28 

PICTURE 

9 

VALUE 

1141 

X  Xt*X 

on A  OA n 
U  U  *T  uo  u 

02 

1ELSW29 

PICTURE 

9 

VALUE 

n 

1  1  Al 

X  X  H  X 

004070 

02 

1ELSW30 

PICTURE 

9 

VALUE 

0. 

1141 

004080 

01 

WORKA. 

1141 

004090 

02 

WORKA-XXXXXGENX. 

1141 

004100 

08 

ASTATNO  PICTURE  9(010). 
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00411C 

08 

1GENDD0013  REDEFINES  ASTATNO. 

1141 

004120 

10 

ASTATNO-XXXXXX0009  PICTURE  9(005). 

1141 

004130 

10 

ASTATN0-XXXXXX0109   PICTURE  9(005). 

1141 

004 140 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0. 

1141 

004150 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0. 

1141 

004160 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0 . 

1141 

0C417C 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+  0. 

1141 

004180 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0. 

1141 

004190 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0. 

1141 

004200 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0 . 

1141 

005000 

08 

FILLER  DISPLAY 

PICTURE   S9( 10)  VALUE 

+0 . 

1141 

005010 

08 

FILLER  DISPLAY 

PICTURE   S9( 10)  VALUE 

+0  . 

1141 

005020 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0 . 

1141 

005030 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0  . 

1141 

005040 

08 

FILLER  DISPLAY 

PICTURE   S9( 10)  VALUE 

+0. 

1141 

005050 

08 

FILLER  DISPLAY 

PICTURE   S9(10)  VALUE 

+0 . 

1141 

00506C 

0  1 

WORKB. 

1141 

005070 

02 

WORKB-XXXXXGENX. 

1141 

005080 

08 

BWATERYR   PICTURE  9(002). 

1141 

005090 

08 

BHIDAY  PICTURE 

9(008) . 

1141 

005100 

08 

BHIDAY3  PICTURE 

9(008) . 

1141 

005110 

08 

BHIDAY7  PICTURE 

9(008) . 

1141 

005120 

08 

BHIDAY15   PICTURE  9(008). 

1141 

005130 

08 

BHIDAY30   PICTURE  9(008). 

1141 

005140 

08 

BHIDAY60  PICTURE  9(008). 

1141 

005 150 

08 

BHIDAY90   PICTURE  9(008). 

1141 

005160 

08 

BHIDAY120   PICTURE  9(007). 

1141 

005170 

08 

BHIDAY183   PICTURE  9(007). 

1141 

005 180 

08 

BHIDAY274  PICTURE  9(007). 

1141 

005190 

08 

FILLER  PICTURE 

X(OOl)    VALUE  SPACE. 

1141 

005200 

01 

IDCLABOOOl. 

1141 

006000 

02 

IDCLABOOOl-GENX 

• 

1141 

006010 

08 

HOLDC  PICTURE  X(OIO)   VALUE  ' 

1141 

00602C 

08 

FILLER  REDEFINES  HOLDC. 

1141 

006030 

14 

HOLD-XXXXXXXXX0009  PICTURE  9(005). 

1141 

006040 

14 

HOLD-XXXXXXXXX0109  PICTURE  9(005). 

1141 

006050 

08 

NINES  PICTURE   9(010)    VALUE  9999999999. 

1141 

006060 

01 

OOIACC. 

1141 

006070 

02 

301ACC  DISPLAY  PICTURE  S9(10). 

1141 

006080 

02 

302ACC  DISPLAY   PICTURE  S9(10). 

1141 

006090 

02 

303ACC  DISPLAY  PICTURE  39(10). 

1141 

006100 

01 

310ACC. 

1141 

006110 

02 

FILLER  PICTURE  X 

• 

1141 

006120 

02 

311ACC. 

1141 

006130 

03 

FILLER  PICTURE  X 

• 

1141 

006140 

03 

312ACC. 

1141 

006150 

04 

FILLER  PICTURE  X 

• 

1141 

006160 

04 

313ACC. 

1141 

006170 

05 

FILLER  PICTURE  X 

• 

1141 

16  of  18 


006180  05  314ACC  PICTURE   X.  1141 

006190  01  320ACC.  1141 

006200  02  FILLER  PICTURE  X.  1141 

007000  02  321ACC.  1141 

007010  03   FILLER  PICTURE   X.  1141 

007020  03   322ACC.  1141 

007030  04  FILLER  PICTURE   X.  1141 

007040  04  323ACC.  1141 

007050  05  FILLER  PICTURE  X.  1141 

007060  05  324ACC  PICTURE   X.  1141 

007070  01  330ACC.  1141 

007080  02  FILLER  PICTURE  X.  1141 

007090  02  331ACC.  1141 

007100  03  FILLER  PICTURE  X.  1141 

007110  03  332ACC.  1141 

007120  04  FILLER  PICTURE  X.  1141 

007130  04  333ACC.  1141 

007140  05  FILLER  PICTURE  X.  1141 

007150  05  334ACC  PICTURE   X.  1141 

007160  01  350ACC.  1141 

007170  02   FILLER  PICTURE  X.  1141 

007180  02  351ACC.  1141 

007190  03  FILLER  PICTURE   X.  1141 

007200  03  352ACC.  1141 

008000  04  FILLER  PICTURE   X.  1141 

008010  04  353ACC.  1141 

008020  05  FILLER  PICTURE   X.  1141 

008030  05  354ACC  PICTURE   X.  1141 

008040  01           STATNO-TEMP,  1141 

008050  02  FILLER  PICTURE   X.  1141 

00806C  02  STATNO-TEMPOIXXXXX  PICTURE  9(005).  1141 

008070  02  STATNO-TEMP02XXXXX  PICTURE  9(005).  1141 

008080  PROCEDURE  DIVISION.  1141 

008090  MOVE  0  TO  lGENERELSW-001 .  1141 

008100  OPEN   INPUT  A114110  OUTPUT  A114151  MESSAGE.  1141 

008110  READX.  1141 

008120  READ  A114110                 AT  END  GO  TO  FINISH.  1141 

008130  IF  lGENERELSW-001  =  1  GO  TO  CKSTAT  ELSE  MOVE  1  TO  1141 

008140  lGENERELSW-001.  1141 

008150  NEWNO.  1141 

008160  MOVE  STATNO  TO  STATNO-TEMP.  1141 

008170  MOVE   STATNO-TEMPOIXXXXX  TO  ASTATNO-XXXXXX0009,  1141 

008180  H0LD-XXXXXXXXX0009.  1141 

008190  MOVE  STATN0-TEMP02XXXXX  TO  ASTATNO-XXXXXXO 109,  1141 

008200  H0LD-XXXXXXXXX0109.  1141 
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C09000  WRITE  OUTPUTX  FROM  WORKA.  1141 

009010  GQ  TO  DETAILS.  1141 

0C9020  CKSTAT.  1141 

009030  MOVE   STATNO  TO  STATNO-TEMP.  1141 

009040  IF   STATNO-TEMPOIXXXXX           =  H0LD-XXXXXXXXX0009  GO  TO  LASTCK.  1141 

009050  IF   STATNO-TEMPOIXXXXX  *    (-1)        =  HOLD-XXXXXXXXX0009  GO  TO  1141 

009060  LASTCK   ELSE  GO  TO   NEWNO.  1141 

009070  LASTCK.  1141 

009080  IF   STATN0-TEMP02XXXXX  NOT  =  H0LD-XXXXXXXXX0109  GO  TO  NEWNO.  1141 

009090  DETAILS.  1141 

009100  MOVE  WATERYR  TO  BWATERYR.  1141 

009110  MOVE   HIDAY  TO  BHIDAY.  1141 

009120  MOVE   HIDAY3  TO  BHIDAY3.  1141 

009130  MOVE  HIDAY7  TO  BHIDAY7.                                                              ■  1141 

009140  MOVE   HIDAY15  TO  BHIDAY15.  1141 

009150  MOVE  HIDAY30  TO  BHIDAY30.  1141 

009160  MOVE   HIDAY60  TO  BHIDAY60.  1141 

00917C  MOVE   HIDAY90  TO  BHIDAY90.  1141 

009180  MOVE   HIDAY120  TO  BHIDAY120.  1141 

009190  MOVE  HIDAY183  TO  BHIDAY183,  1141 

009200  MOVE   HIDAY274  TO  BHI0AY274.  1141 

OlOOOC  WRITE   OUTPUTX  FROM  WORKB.  1141 

010010  GO  TO  READX.  1141 

010020  FINISH.  1141 

010030  MOVE   NINES  TO  ASTATNO.  1141 

010040  WRITE   OUTPUTX  FROM  WORKA.  1141 

010050  CLOSE  A114110  A114151.  1141 

C10060  MOVE    • 1114101N0RMAL   END  OF   JOB'    TO  MES-TM.  1141 

010070  WRITE   MES-TM.  1141 

010080  CLOSE  MESSAGE.  1141 

010090  STOP  RUN.  1141 


C  FORMATS  DATA  FROM  GS 

C  JOB  NUMBER  920603-0030 

C  PROGRAM  NUMBER  1183 

DIMENSION  X(10) 
1   READ(1,6)  K,J,N 
6  FORMAT( 12, I8,1X,I8) 
IF   (K-99)  32,30»30 
30   IF   (J-99999999)  42t20,20 
32  IF   (K-l)  34,42,42 
34  IF   (N-1)  36,42,42 
36  WRITE(6,38)  K,J 
38  FORMATdHl,  12,18) 

WRITE(4,40)  K,J 
40  F0RMAT(1HV,I2,I8) 
GO  TO  1 
42         BACKSPACE  1 

2  READ(    1,4)    K,(X(I ),I=1,10) 

4  FORMAT ( I2,7F8.2,3F7.2) 

WRITE(6,16)  K,(X(I),I=1,10) 
16  FORMAT( IHO, I2,7F8.1,3F7. 1) 

WRITE(4,54)   K,(X(I ) ,1=1,10) 
54  FORMATdHV,  I2,7F8.1,3F7.  1) 

GO  TO  1 
2D  REWIND  1 

WRITE  (6,99999) 
99999  FORMAT( 1H1,20X,10HEND  OF  JOB/lHl) 
STOP  999 
END 


